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Abstract 

SUSY models with a gauge singlet easily allow for a strong first order electroweak 
phase transition (EWPT) if the vevs of the singlet and Higgs fields are of comparable 
size. We discuss the profile of the stationary expanding bubble wall and CP-violation 
in the effective potential, in particular transitional CP-violation inside the bubble 
wall during the EWPT. The dispersion relations for charginos contain CP-violating 
terms in the WKB approximation. These enter as source terms in the Boltzmann 
equations for the (particle-antiparticle) chemical potentials and fuel the creation of 
a baryon asymmetry through the weak sphaleron in the hot phase. This is worked 
out for concrete parameters. 
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1 Introduction 



Different from models starting at the Grand Unified scale, the ingredients for a creation 
of the baryon asymmetry of the universe in a first order electroweak phase transition 
(PT) have a fair chance to be tested experimentally in the near future. Besides non- 
equilibrium the other two necessary criteria of Sakharov for baryogenesis, CP and baryon 
number violation are naturally fulfilled by the CP violating bubble wall of a first order PT. 
The bubble wall separates the hot symmetric phase with rapid sphaleron transitions from 
the Higgs phase, where sphaleron transitions are suppressed by exp(v(t)/T). However it 
turned out that the electroweak Standard Model (SM) does not have a strong first order 
PT, indeed there is a crossover behavior and no PT at all for Higgs masses beyond the 
present experimental bound [I]. Furthermore, CP- violation by the CKM matrix is very 
small. In spite of tremendous successes of the SM it is commonly believed that it has 
to be embedded / enlarged to a more fundamental theory. Supersymmetry is supposed 
to be an important facet of such a theory but it is still not backed by experiments. In 
concrete SUSY models, in order to have a strong first order PT one has to strengthen 
the (one loop) u tp 3 " term of the effective potential or to have such a term already on 
the tree level. In the minimal SUSY extension of the SM, the MSSM, the superpartner 
of the right-handed top - with a mass below that of the top - gives such a strong (p 3 
loop correction. With a negative SUSY breaking scalar mass it is almost massless in the 
hot phase and has a strong Yukawa coupling to the Higgs fields. This is confirmed by 
detailed analytical and lattice calculations ||, showing a strong first order PT even for 
Higgs masses as large as 110 GeV just beyond the experimental bound for MSSM Higgses. 
There are also possible CP-violating phases in the Higgs Lagrangian if 1-loop effects are 
taken into account 0. They are restricted much by the measurements of the neutron 
electric dipole moment unless one invokes special conditions. A recent investigation of a 
possible spontaneous CP- violation just in the bubble wall at the temperatures of the PT 
gave negative results 0. 

In this paperf] we want to discuss baryogenesis in a strong first order electroweak PT 
(EWPT) in much detail starting from a model where all the three criteria given above can 
be fulfilled without much problem. This is a supersymmetric model with an additional 
gauge singlet superfield. In its original form (NMSSM) is was designed to substitute the 
problematic \iH\H 2 term of the MSSM superpotential by a coupling SH\H 2 @ 

W Z3 = \SHtH2 + ^S 3 

with a ^-symmetry whose spontaneous breaking causes dangerous domain walls in the 
early cosmos [0]. The soft SUSY breaking potential contains a term \A\SHiH 2 which 
can act as a if 3 term if both the vevs of the Higgses and of the singlet are of the same 
order of magnitude [H. But in a Z 3 -symmetric model with universal SUSY breaking this 
turns out not to be possible if one wants to have a reasonable spectrum of particles. 

In agreement with ref. || we thus add further terms to the superpotential which now 
obtains the form 

W = XSH 1 H 2 + ^S 3 + fiHiHz + rS (1.1) 

3 whose content was presented at the TMR network meeting Finite Temperature Phase Transitions in 
Particle Physics in Korfu, September 1999. 
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and is not Z 3 -symmetric anymore. The soft SUSY breaking Lagrangian contains scalar 
masses, gaugino masses and A-terms 



A oft = \A x SHiH 2 + ^A k S 3 + Y e A e e c lH 1 + Y d A d d c ~qH x + Y u A u u c qH 2 + h.c. (1.2) 

Thus we have reintroduced a /x-term and we have the associated fine-tuning problem 
known from the MSSM l lQfl . Additionally, there is the danger of quadratically divergent 



singlet tadpoles [|TT],|T^]. Such tadpole diagrams require three ingredients: (i) a singlet 
field (ii) non-renormalizable interactions and (iii) soft SUSY breaking terms, (i) Giv- 
ing the "singlet" a charge under some discrete symmetry one can remove the tadpoles 
from the very beginning. But the discrete symmetry cannot be exact in order to avoid 
unacceptable domain walls at the EWPT [13]. (ii) In refs. [12,14] models with gauged 



R-symmetry or duality symmetry, both broken at some superheavy scale have been pro- 
posed to forbid dangerous non-renormalizable operators. On the renormalizable level 
these models typically have no discrete symmetries, and therefore are not plagued by the 
domain wall problem. While such models do not solve the /z-problem, they are still very 
interesting with respect to Higgs phenomenology and the EWPT. For that reason we will 
study this type of singlet model in the following, (iii) Another way to evade the tadpole 
divergences restricts the soft SUSY breaking terms. "Gauge mediated SUSY breaking" 
(GMSB) in the context of singlet models does not have domain wall problems. A \i- 
parameter is generated by radiative corrections and the singlet vev |15|]. However one of 
the special properties of GMSB models seems to be the strong suppression of A-terms. 
These however also contain the </2-type term responsible for a strong first order PT. 

In order not to have too many parameters we use universality of SUSY breaking at the 
GUT scale and run renormalization group equations to get down to the electroweak scale. 



With such a Lagrangian and adding the temperature corrections we can demonstrate 
that one can easily get a strong first order PT for Higgs masses as high as 115 GeV and 
superpartner masses beyond present experimental limits. The effective potential in H\, 
Hi and S can be used to derive differential equations for the bubble wall profile and for 
varying CP- violating phases in the bubble wall arising from constant explicit phases in the 
theory or, more interesting, from spontaneous CP-violation. This is described in sections 
4 and 5. As worked out in section 6, in the WKB approximation the CP-violating phases 
in the bubble wall create terms in the dispersion relations which differ between particles 
and anti-particles. This gives rise to a driving term in the Boltzmann equation for the 
difference of particle and anti-particle chemical potentials, in particular for left-handed 
quarks and their superpartners. This difference is converted into a baryon asymmetry by 
the hot phase sphalerons in front of the proceeding bubbles. 



2 The model 

Weak scale supersymmetric models have many unknown parameters related to super- 
symmetry breaking. A considerable part of this parameter space is excluded because of 
FCNC or the appearance of charge and color breaking vacua. These problems are partially 
evaded by imposing universality of the soft terms at the GUT scale, Mqut ~ 2.6 • 10 16 
GeV, where the SM gauge couplings unify. At the GUT scale there is a common gaugino 
mass M , a universal scalar mass squared Wq and a universal trilinear coupling A . Thus, 



2 



universal SUSY breaking drastically reduces the parameters of the model to 

yto, A , fco, M , Aq, mjj, /i , r , B . (2.3) 

The subscript "0" indicates that the masses and couplings are evaluated at the GUT 
scale. y t denotes the top Yukawa couplingfj, and \iBH x H 2 is the soft Higgs mass term 
corresponding to the //-term in the superpotential. Since the masses of the top quark and 
the Z-boson mass are known, "only" seven parameters of ( |2.3| ) are independent. 

To evolve the parameters from the GUT scale to the weak scale we use the renormal- 
ization group equations (RGEs) in the 1-loop approximation. For the ^-symmetric case 



the relevant RGEs have already been given in ref. [17]. Defining t = In Q the ^-breaking 



terms /x, r and B are easily included with help of ref. [JIS 

j t » = ^^y 2 t +2X 2 + 2k 2 -?>gl-gl) 



d 1 
It ~ 16^ 2 ' 



d „ 1 



il B = T^(^ 2 B + 3y 2 A t + 2X 2 A x + 3g 2 2 M 2 + g 2 1 M 1 ). (2.4) 
dt 1d7t z 

Even with the simple universal pattern of soft breaking terms we are left with a 9- 
dimensional parameter space. Additionally, we must satisfy two constraints coming from 
the Z-boson and top quark masses. In the literature (see e.g. ref. [|HJ), usually "random 
shooting" is applied to deal with such a situation. This means randomly chosen sets 
of GUT scale parameters are evolved to the electroweak scale, where their phenomeno- 
logical implications are investigated. Although the physical Z-boson mass can easily be 
reproduced by an appropriate rescaling of the dimensionful parameters, one typically is 
still plagued by some light unobserved SUSY particles or an unphysical top quark mass. 
Therefore this procedure is rather inefficient. Furthermore, random shooting only pro- 
vides statistical averages and correlations. To avoid these shortcomings we use the more 
systematic approach which we introduced in ref. ||16|| . It allows the elimination of /i, r and 



B by the Higgs and singlet vevs, while maintaining universal SUSY breaking. It relies 
on the observation that the ^-breaking terms /i, r and B do not enter the RGEs for the 
remaining parameters. They can therefore be calculated without specifying the former 
ones. Our procedure consists of the following steps and is also sketched in fig. [[]: 

• Fix the values of x = (S), tan/5 = v 2 /v\, A and k at the weak scale. The choice of 
tan/3 automatically fixes y t by the relation y t vsm/3 = m t0 p{= 175 GeV). 

• Evolve y t , A and k to the GUT scale. 

• Choose a set of the parameters Mq, Aq and m\ at the GUT scale. 

• Run all parameters down to the weak scale, with exception of /i, r and B whose 
initial conditions have not yet been specified at this point. 

• At the weak scale calculate /i, r and B by using the saddle point conditions for 
the 1-loop Higgs potential d I V(H 1 ,H 2 ,S) = (I = H u H 2 , S) at tff = H+ = 0, 
#1,2 = u i,2, S = x. 



4 We neglect the Yukawa couplings of the leptons and light quarks. This approximation is well justified 
in the range of small and mediate tan (3 which we consider in the following. 
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Figure 1: Sketch of our procedure to fix the weak scale parameters using RGEs and saddle 
point conditions. 



Finally one can use the corresponding RGEs to determine the values of //, r and B at 
the GUT scale. Appropriate values of fi , r and B can always be found because of the 
effectively linear structure of eqs. (2^). Our procedure has the additional benefit that the 
constraints coming from the top quark and Z-boson masses, M§ = (g\ + g%){v\ + vD/Z, 
are automatically built in. This reduces the number of free parameters by two so we are 
finally left with the 7-dimensional parameter space 



tan p, x, A, k, M , A , m . 



(2.5) 



Of course, not every parameter set leads to a viable phenomenology. The constraints 
on the model parameters are related to the particle spectrum and the vacuum structure. In 
the elimination procedure discussed before we assumed that the extremum parametrized 
by v, tan/? and x is indeed the global minimum of the scalar potential. Thus we must 
(numerically) verify that there are no deeper minima in V(Hi, H2, S). If there appears a 
deeper minimum the parameter set has to be discarded. 

Additional constraints arise from the required absence of charge and color breaking 
minima (CCB minima) deeper than the standard minimum, i.e. from the absence of squark 
and slepton vevs. We check slepton and squark vevs induced by large trilinear couplings 
|2(||21]]. Furthermore, there are CCB minima which come from negative scalar mass 
squares, typically m% . These dangerous directions involve Higgs, squark and slepton 
fields ("UFB" directions). The condition for such a minimum not to be deeper than 
the standard minimum implies a lower limit on the ratio mg/Mp of 0(1) ||22|| . However, 
the decay rate of the standard vacuum into a minimum in the UFB direction is usually 
negligible compared to the age of the universe. We therefore allow for a meta-stable 
standard vacuum with respect to the UFB directions and disregard the corresponding 
constraints. 

Up to now, no superpartners of the SM particles have been detected. From the experi- 
mental lower limits on the SUSY particle masses |23j various constraints on the parameter 
space can be derived. The experimental limit on the chargino mass M~± > 90 GeV 

Xi — 

translates into bounds on the universal gaugino mass and the /1 parameter 



I M 1 > 100 GeV, |/i + Xx\ > 80 GeV. 



(2.6) 
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Figure 2: Scan of the M -A plane for two different sets of (x, tan/3, k). The colored 
regions represent the phenomenologically viable range of the parameters before cuts from 
the Higgs boson search are applied. The dashed lines are curves of constant mass of 
the lightest CP-even Higgs boson. In the blue (green) areas, the lightest Higgs boson is 
predominantly a singlet (Higgs) state. 



If the bound on Mo is satisfied the gluino mass M3 is automatically above its experimental 
limit f23[. In the case of mo > 100 GeV also the squark and slepton masses are compatible 
with the experimental data. 

Of particular importance are the properties of the lightest neutral CP-even 'Higgs' 
mass eigenstate h, which is a mixture of the Higgses and the singlet. If h has a large 
singlet content its coupling to the Z-boson and thus its production cross section at LEP is 
significantly reduced . Parameter sets ( |2.5| ) cannot be ruled out by simply calculating 
the lightest Higgs mass M h . For example, if the Higgs production cross section is reduced 
by a factor of 10 compared to the SM, Higgs masses down to about 60 GeV are still in 
agreement with the data ||26|| . 

In fig. § we present a scan of the Mq-Aq plane for two different sets of (x, tan /3, k) fllfl . 
The parameters M and A are of particular interest since they determine the values of 
the trilinear couplings at the weak scale which will play a prominent role in the discussion 
of the EWPT in section 3. We have taken m > 100 GeV in order avoid light squarks 
and sleptons. Choosing a small value of A and a large value of k prevents the appearance 
of very light Higgs bosons. In fig. |2|b we increased tan/3 and the singlet vev to obtain 
large Higgs masses up to 115 GeV. The increase of M^ with decreasing A can be traced 
to the singlet diagonal entry in the Higgs boson mass matrix which is diminished by the 
Ak contribution. On the same lines the increasing singlet content of the lightest CP-even 
Higgs state with increasing A can be understood. In fig. 0b this effect is reduced by the 
larger values of \x\ and k, which render the singlet state rather heavy. 

The most important constraints on M and A arise from the chargino mass and 
from the vacuum structure. The lower bound on M is a consequence of eq. ( |2.6| ). The 
requirement of the standard minimum being the global minimum of the Higgs potential 
implies the upper and lower bounds on Aq displayed in fig. |2|. Constraints from the 
required absence of A-term induced CCB minima are then automatically satisfied. Note 
that M is also bounded from above. In fig. || we excluded parameter sets which predict 
M h < 65 GeV. This leads to the upper bound on A in fig. 0a. In fig. no additional 
constraint arises, because of the larger Higgs masses. On the other hand, if we would 
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allow for smaller Higgs masses in fig. [|a, still an upper bound on A Q would be implied 
by the vacuum structure, as is the case in fig. |2|b. Taking into account the data from 
the Higgs boson search further restricts the upper bounds on M and A . Roughly, the 
regions M > 600 GeV in the tan/3 = -5 set, and M > 1200 GeV in the tan/3 = -10 
set are additionally excluded by the Higgs boson data ||26|| . These parameter regions are 
characterized by a very heavy spectrum of SUSY particles. For that reason they are not 
very promising candidates for the implementation of electroweak baryogenesis anyway, as 
will be discussed in section 6. 

We close this section by emphasizing some important differences between our model 
and the ^-symmetric NMSSM, which is usually considered in the literature, e.g. in 



refs. [[L9|,[25|]. By setting fj, — eq. ( |2.6| ) can be translated into a lower bound on the 
singlet vev, \x\ > 125 GeV. However, the analysis of the ^-symmetric case carried out in 
ref. [[B|] shows that the actual lower bound on \x\ is much larger, once all phenomenolog- 
ical constraints are taken into account. Typically one has a singlet vev in the multi-TeV 
range and couplings A < 1 and k <C 1, with Ax and kx = 0{Mz). In the next section 
we will discuss that in such a scenario the singlet is just a spectator during the EWPT, 
which proceeds more or less in the standard (i.e. MSSM) way. This property additionally 
motivates our inclusion of the //-term in the model, while it seems questionable whether 
the idea of electroweak baryogenesis can be realized in the Z 3 -symmetric NMSSM. An 
additional peculiarity of the ^-symmetric model is that large A-terms are required for 
successful electroweak symmetry breaking, leading to A\ > 9mQ. However, this param- 



eter range is severely constrained from the absence of CCB minima [21 1. Some of the 



special features reported in ref. [19] may be due to the assumption of universal soft SUSY 



breaking used in that work. It seems questionable, however, if the problematic bound, 
^> Mz, can be evaded in the case of non- universal soft terms [pTfl , since the required 



large values of A typically lead to small values of the Higgs mass due to large Higgs singlet 
mixing. The NMSSM with and without Z 3 symmetry also have different properties with 
respect to CP-violation, as will be discussed in section 5. 



3 Strength of the electroweak phase transition 

Order and strength of the EWPT are central questions in electroweak baryogenesis. Only 
in the case of a first order phase transition (PT) the associated departure from equilibrium 
is sufficient to induce a relevant baryon number production. To avoid baryon number 
washout after the PT the even stronger criterion v c /T c > 1 has to be satisfied, where v c 
denotes the Higgs vev at the critical temperature T c . A first order phase transition is 
triggered by cubic terms in the finite temperature effective potential [^8 |. In the (MS)SM 
these terms arise from 1-loop thermal corrections of bosons and therefore are small from 
the very beginning. Thus it is difficult to satisfy v c /T c > 1. In the NMSSM, on the other 
hand, trilinear terms enter already the tree-level Higgs potential due to Higgs singlet 
couplings, leading to a significantly stronger EWPT |8|, p|, |I6||. These contributions stem 
from the soft SUSY breaking couplings Ax and A^, and from the //-term 

(\fi*S + h.c.)(|#iT + |#°| 2 ) + (XAxSH ^ + ^A k S 3 + h.c). (3.7) 

All these trilinear terms explicitly contain the singlet field. Therefore, in order to induce 
deviations from the (MS)SM behavior also the singlet vev must change during the phase 
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transition. Since at the EWPT thermal contributions to the effective potential are of the 
order 0(T 4 ) ~ M|, this requires the mass and the vev of the singlet also to be of the 
order of the electroweak scale. 

At finite temperature the effective potential of the neutral Higgs and singlet fields is 
modified by the interaction with the hot plasma 

V T (H°, Hi S) = V trcc {Hl Hi S) + VW(HI H° 2 , S) + v£\hI H° 2 , S). (3.8) 

In the 1-loop zero temperature corrections we include tops, stops and gauge bosons, 
while in the 1-loop finite temperature part V T ^ also Higgs bosons, neutralinos and charginos 
are taken into account. We do not make a high temperature expansion, as some of the 
particles can be heavy in part of the field space. Rather, V^ 1 -* is evaluated using a spline 
interpolation between the high and low temperature regions. 

We stress that in the range of parameters we will study a strong first order EWPT is 
the consequence of the tree-level terms in the Higgs potential fl3.7p . Since the stop mass 
will turn out to be always larger than 200 GeV, the thermally induced cubic terms are 
too small to account for v c /T c > 1 0. The most important finite temperature effect is 
the appearance of thermal effective masses 

m 2 — > m 2 (T) = m 2 (T = 0) + const • T 2 , (3.9) 

where the constant encodes the couplings of the Higgs and singlet fields to the particles 
in the plasma. It is this positive thermal contribution to the Higgs mass that makes 
the symmetric phase stable at high temperatures, and thus causes the restoration of the 
electroweak symmetry. Using the complete 1-loop expression rather than the simple high 
temperature approximation, is just a convenient prescription to handle the decoupling of 
heavy particles. Because of the dominance of the tree potential, we neglect contributions 
to the thermal potential stemming from daisy resummation |29],|30| and 2- loop diagrams 
I- 

In order to determine the strength of the PT one has to compute the critical temper- 
ature T c . We define T c as being the temperature where the symmetric and the broken 
minimum of the finite temperature Higgs potential (|3.8|) become degenerate. We denote 
the Higgs and singlet vevs in the broken phase by (v i c , t>2 C , x c ). In general, the singlet 
vev is different from zero even in the symmetric phase. We refer to a PT as being 
"strongly first order", if it avoids baryon number washout, according to the condition 
v c /T c > 1. Since both Higgs fields contribute to the sphaleron energy, v c is given by 

v c = a/2 (|^ic| 2 + |^i c | 2 ) 1//2 , where the factor of \pl is due to our choice of field normaliza- 
tion. Although CP-violation is essential for baryogenesis, we assume CP-conservation in 
this section, i.e. all vevs, mass parameters and coupling constants are taken real valued. 
Of course, we have to verify that Vr has no deeper CP-violating minimum. Turning on 
phases much smaller than one would induce only marginal changes in our results. 

In this paragraph we study again the two parameter sets already discussed in the con- 
text of fig. H in section 2. We determine the critical temperature by numerical minimiza- 
tion of Vt- After having checked that no CP- violating minima are present we disregard 
the imaginary parts of the Higgs and singlet fields.^ We then only have to minimize Vr 
in the real valued fields Re(if°), Re(H 2 ) an d Re(S'). 

5 Bubble walls in the presence of CP-violation will be discussed in section 5. 
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Figure 3: Scan of the M Q -A plane for the two different sets of (x, tan/3, k) which have 
already been considered in fig. |2|. In the red (yellow) areas the PT is strongly (weakly) 
first order, i.e. v c /T c > 1 (v c /T c < 1). The dotted lines are curves of constant mass of the 
lightest CP-even Higgs boson. In the regions below (above) the dashed lines the lightest 
Higgs boson is predominantly a Higgs (singlet) state. 



Our investigations of the strength of the PT are summarized in fig. |^ where the regions 
of strong and weak PT in the Mq-Aq plane are displayed for the two parameter sets of 
fig. |2| fL6 |. In case of tan/3 = —5 and x = —100 GeV (fig. |3]a) the PT is strongly first 



order in most part of the phenomenologically allowed range of parameters. However, the 
corresponding Higgs masses up to 90 GeV are compatible with the experimental Higgs 
mass bounds only because of the reduced Higgs production cross section due to Higgs 
singlet mixing. In the parameter set of fig. ^|b we increased | tan /5| and \x\ in order to 
obtain larger Higgs masses. As a consequence, the region of weakly first order PT is 
enlarged. However, a strong PT occurs still for a wide range of the parameters, while 
Higgs masses up to 115 GeV are consistent with v c /T c > 1. 

In contrast to the SM or the MSSM where the PT definitely becomes weaker with 
increasing Higgs masses the situation in the NMSSM is more complicated. Larger Higgs 
masses can be related to a stronger phase transition. This may happen for example if the 
broken and symmetric minima of the effective potential are almost degenerate already at 
zero temperature. The strong PT at negative A and M h > 100 GeV in fig. |3|b results 
precisely from this effect. 

Our results |16| are in reasonable agreement with those of ref. 0, where the strength 
of the EWPT in the general NMSSM has also been studied. The authors found that 
about 50% of their parameter sets were compatible with v c /T c > 1, and also advocated 
for x ~ v in order to obtain a strong phase transition. While in ref. |J no restrictions with 
respect to the SUSY breaking were made, we demonstrated that even with universal SUSY 
breaking a strongly first order PT is quite natural in the general NMSSM. Furthermore, we 
used updated (more restrictive) experimental bounds to constrain the parameter space. 
A further difference comes from the procedure used in scanning the parameter space. 
In ref. |5J random shooting was applied, which only allowed for statistical statements. 
Moreover, from the initially tested 105000 parameter sets just 2% gave a viable zero 
temperature phenomenology, leaving merely about 2000 sets for the study of the PT. This 
again demonstrates the usefulness of the systematic procedure for scanning the parameter 
space, which we described in section 2. 
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4 Shape of the phase boundary 



In the previous section we defined T c as being the temperature where the symmetric and 
the broken minimum of Vr become degenerate. However, tunneling with formation of 
bubbles of the broken phase starts at some lower "nucleation" temperature, when the 
symmetric phase is already meta-stable. The bubble wall profile will crucially enter the 
calculation of the baryon production during the PT which will be discussed in section 
6. Here we concentrate on CP-conserving bubble walls. The very important case of 
CP-violating wall profiles will be discussed in the next section. 

At high temperatures the probability for thermal tunneling is proportional to e~ 53//T , 
where S3 is the three-dimensional action of the static field configuration describing 



tunneling pi||. Here $ collectively denotes the Higgs and singlet fields. Assuming spherical 



symmetry the bubble configuration ( "critical bubble" ) obeys the equation of motion 

dr 2 r dr <9$ 

Further simplifications of eq. ( }4.10|) are justified if the tunneling occurs between two 
almost degenerate minima of the potential with an energy difference AVr much smaller 
than the height of the potential barrier. In such a case the radius of the bubble becomes 
much larger than the thickness of the bubble wall, which thus is referred to as "thin wall 
limit". Neglecting therefore the dQ/dr term in ( |4.10| ) we are left with 



d 2 <$> d 

^-8,VH»)=0. (4.11) 



We have replaced the spatial coordinate r by z, indicating that the solution to ([4.1 1|) 



may be viewed as a planar domain wall with translational invariance in the directions 
perpendicular to the z axis. During the period of stationary expansion the pressure 
induced by the energy difference of the minima, AVr, is compensated by friction |67|,|68[ . 



In the following we model this effect by taking the effective potential at the critical 
temperature. 

Notice that ( |4 . 1 1| ) is just the classical equation of motion of a particle moving from 
one maximum of the turned around potential — V(Q) to the other, where is comes at rest. 
In this picture z takes the role of the time coordinate and $ represents the configuration 
space variable. Obviously, this a very delicate process, especially if more than one scalar 
field is involved. Small changes of the initial conditions lead to a completely different 
shape of the trajectory. 

For a general effective potential the bubble wall equations have to be solved numer- 
ically. In the case of only one scalar field the so called "overshooting undershooting 
procedure" can by applied: The initial value $(r ) is tuned such that the trajectory 
approaches $ = •I'sym in the limit r — > 00, which then gives the bubble shape. This 
procedure can be used for the critical bubble as well as for the domain wall configuration. 

The situation is completely different once there are additional directions in field space. 
Although in principle the shooting procedure is still applicable, in practice the required 
high accuracy of the initial conditions cannot be achieved. Thus one has to devise other 
numerical methods P,|33|-|37|. They determine the bubble wall configuration as the min- 



imum of a functional J-"[&] which is built from the squared equations of motion [35] 



Details concerning the numerical algorithms used to minimize J 7 on a lattice are given in 
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refs. |33|,|35|]. It turns out that the minimization procedure is reliable only if the starting 
configuration, i.e. the ansatz for the wall shape, differs not too much from the actual 
solution. Otherwise one keeps being stuck to local minima, which arise from discretizing 



the functionals on a lattice [Q. Thus minimization starting from an arbitrary initial con- 
figuration is not possible up to now, and finding an appropriate ansatz is very important. 
In the case of the domain wall the kink solution fits the exact bubble shape reasonably 
well, especially if every field is allowed for having its own wall thickness 

0^) = |(l-tanh(^)V (4.12) 



In general, different off-sets Si which shift the fields against each other are possible as 
well. For the domain wall case the mechanical analogue to energy conservation, E(z) = 
~(d$/dz) 2 — V(<&) = const, provides a very sensitive criterion to check the quality of the 
numerically obtained solution. 

In general, the numerical methods developed in refs. [|5],[33-36j are inevitable to de- 



termine the bubble shape in multi-field models. But if the problem is effectively one 
dimensional some "improved shooting method" can still be applied. The investigation of 
CP-conserving bubble wall shapes in the MSSM [0,^] revealed that for realistic Higgs 
masses the variation of tan/3 in the bubble is very small, 5(3 ~ 0(1O~~ 2 — 1CT 3 ). Thus the 
bubble is very accurately described by taking only the combination of the Higgs fields, 
H, which corresponds to the direction of the broken minimum. Simple shooting along 
this direction provides excellent values for the wall thickness, the surface tension or the 
action of the critical bubble. Even the small variation in f3 can be reliably determined 
by minimizing the effective potential along the direction perpendicular to H . The line 
which connects these minima gives a good approximation to the trajectory in field space 
the actual bubble solution is corresponding to. This can be checked by using the "energy 
conservation" criterion. The kink ansatz turns out to be a good fit for bubble profile. 

In case of the NMSSM the straight connection between the broken and the symmet- 
ric minimum is no longer a reasonable approximation to the true bubble wall trajectory. 
Ignoring the variation of tan j3, which turns out to be small and hence can be included af- 
terwards, we effectively have a two field system, H and S. In a large part of the parameter 
space (as long as A Q is significantly above its lower bound shown in fig. |2|) the effective po- 
tential is characterized by a distinct smooth ridge which connects the symmetric and the 
broken minimum. This ridge S = F(H) can be determined by minimizing the potential 
along the direction perpendicular to the straight connection of the symmetric and broken 
minima.^ Using Vr(H, F(H)) the set of bubble wall equations ( |4.11| ) reduces to only one 



differential equation which can again be solved via the shooting method. However, not 
every function F is a possible candidate for the actual bubble wall. Since dH/dz vanishes 
in the symmetric and broken phase, integrating the bubble wall equation results in the 
constraint 



dH—V T (H,F(H))=0. (4.13) 



6 Numerically, the following method leads to better results: Determine the saddle point along the ridge, 
i.e. BhVt = dsVr = 0. The ridge is then defined as the trajectory of an over-damped particle rolling 
from the saddle point down to the (symmetric or broken) minimum, satisfying d$/dz + d$Vr = 0. 
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Figure 4: (a): Bubble wall profile at the critical temperature for the parameter set of 
fig. 0a with M = 300 GeV and A = 0. (b): The corresponding trajectory in the H- 
S plane (solid line) and the straight connection between the symmetric and the broken 
minimum (dashed line). (All units in GeV.) 



In order to evade this inconsistency we deform our ansatz by F(H)=F(H)+fH(H-v c ), 
where the free parameter / is fixed by eq. ( [4.13j ). This ansatz is motivated by the me- 
chanical analogue of the bubble wall problem: if a ball rolls from one maximum of — 
to the other, it moves a bit below the ridge in order to compensate for the centrifugal 
force. On half way between the two maxima the centrifugal force is most efficient and the 
trajectory F deviates maximally from the ridge F. In practice, F — F is small compared 
to the deviation of F from the straight line. 

Instead of expressing S in terms of F(H) one could just as well eliminate H by F~ 1 (S). 
If F corresponds to the trajectory of the actual bubble configuration the two prescriptions 
are equivalent. For a general F the two prescriptions will lead to different results for the 
wall profile. The difference between the two solution indicates the quality of the chosen 
F. We emphasize that although several steps are necessary to compute the bubble wall 
shape via this improved shooting approach, the required numerics is fairly simple. It 
can by carried out using computer algebra systems like Maple or Mathematica. The 
results are competitive to those obtained with the sophisticated minimization algorithms 
of refs. |33|,|35|], which on the other hand were very helpful to check the validity of our 



approximations. 

Let us discuss a specific example from the parameter set of fig. 0a with M = 300 
GeV and A = 0. It corresponds to a lightest Higgs mass = 80 GeV, which because 
of the reduced coupling to the Z-boson is compatible with the experimental data. We 
have T c = 109 GeV and v c /T c = 1.12, i.e. the washout of baryon number after the phase 
transition is avoided. Using our improved shooting method we find the bubble wall profile 
displayed in fig. [|a. In fig. [|b we show the corresponding trajectory in field space which 
considerable deviates from a straight line. Fitting the numerical solutions by the kink 
ansatz Q4.12|) we obtain the wall thicknesses for the Higgs and singlet fields, = 0.13 
GeV -1 = 14/T C and L s = 0.10 GeV -1 = 11/T C , respectively. For the surface tension 
we find a = 46300 GeV 3 . The energy conservation check gives AE = 6000 GeV 4 which 
has to be compared to the barrier height V& ~ 125000 GeV 4 . This is the same level of 
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Figure 5: (a): Variation of tan/3 in the bubble wall at the critical temperature for the 
parameter set of fig. [|. (b): Shape of the critical bubble for the same parameter set at 
T = 108.5 GeV. (Units in GeV.) 



accuracy [[37] which is accessible in the minimization approach of ref. [|33|]. On the other 
hand, taking the simplest ansatz for F, the dashed straight line in fig. |]b, results in 
L h = L s = 8/T c and a = 82200 GeV 3 , which is already off by about a factor of two. 

Up to now we neglected the variation of tan (3 in the bubble wall which has a strong 
impact on some sources for baryogenesis, as will be discussed in section 6. In order to 
include this effect in our calculation we simply minimize the effective potential in the 
direction orthogonal to H while keeping H and S fixed, implicitly assuming that it is 
only a small perturbation to the trajectory F(H). In fig. [5]a we show the variation of 
tan/? for the bubble wall of fig. f|. We obtain 5(3 = 1.2 x 10~ 3 , i.e. the assumption of 5(3 
being a small perturbation is very well justified. Our result is in complete agreement p7[] 
with the one computed by using the minimization technique of ref. [Q. Moreover, we 
find that the variation of tan (3 in the NMSSM and MSSM RSfl are of the same order of 
magnitude. Particularly, the singlet field provides no additional sources for 5(3 which can 
be traced to its equal couplings to both Higgs fields. 

Finally, our improved shooting method can be used to determine the shape of the 
critical bubble ( [4.10| ). The first bubbles nucleate when the action of the critical bubble 
satisfies S3 (T n ) /T n ~ 130 — 140 . In the computation of the T n the small tan (3 variation 



in the wall can be safely ignored. Considering again the parameter set of fig. |] we obtain 
the bubble configuration shown in fig. |5|b, where T = 108.5 GeV. Since T c = 109.2 GeV 
the super-cooling amounts to 0.7 GeV, the same order of magnitude as in the MSSM 



35|. Notice that the wall shape is still very similar to a kink solution (|4.12|), which is cut 



off at r = 0. Calculating the corresponding action we get S 3 /T = 134. Thus we are at 
the nucleation temperature. Note that the thin wall approximation 5*3 = i^p<7 3 /(AVr) 2 , 
leading to S3/T = 94 is not very accurate in this case. Here AVr = 4.04 x 10 5 GeV 4 
is the potential barrier. Since in our example bubble radius and wall thickness are of 
comparable size this difference is not surprising. 

We emphasize that our method to determine bubble shapes is restricted to cases where 
the bubble trajectory is associated with a smooth ridge in the effective potential, which 
correspond to values of Aq significantly above its lower bound. Otherwise one is referred 
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to the sophisticated minimization methods of refs. [33, 35,36 . 



5 CP-violating bubble walls 

Having established a strongly first order EWPT in the NMSSM, successful electroweak 
baryogenesis still depends crucially on the available amount of CP-violation. In models 
containing two Higgs doublets, such as the MSSM or the NMSSM, complex vevs of the 
Higgs (and singlet) fields provide additional sources of CP-violation. Complex vevs can 
either arise spontaneously or are induced by explicitly CP-violating couplings. In this 
section we discuss their impact on bubble wall profiles in the NMSSM. Most interestingly, 
we find that CP-violation may be restricted to the phase transition itself, a phenomenon 
which is called "transitional CP-violation" . In particular, this scenario allows for a large 
amount of CP-violation during the process of baryon production, while being completely 
unconstrained by experiment. 

To start with, let us summarize the case of CP-violation at zero temperature. In the 
MSSM the tree-level Higgs potential automatically conserves CP. However, explicit CP- 
violation emerges from possible complex phases of the soft SUSY breaking A-terms and 
gaugino masses, and from the /z-parameter Q. These phases contribute to the electric 



dipole moments (EDMs) of quarks and electrons |38] . From the experimental upper limits 



HI on the neutron EDM, d n < 1.1 x l(T 25 e cm, and the electron EDM, d e < 4.3 x l(T 27 e 



cm, constraints on the supersymmetric phases can be derived. Conservatively, rather small 
super symmetric phases of the order O(10~ 2 — 1CT 3 ) are required to satisfy the experimental 
bounds P^] . Larger phases may only be tolerated if the first and second generation squarks 



have masses in the TeV range ||40|| , or if accidental cancellations do occur. Recently it was 



realized that these cancellations are more generic than thought previously ||41|| . Explicit 



CP-violation occurring in the SUSY breaking sector of the MSSM is communicated to the 
Higgs sector by radiative corrections. However, the complex phases induced in the Higgs 
vevs are much too small to have any phenomenological implications ||42|| . Nevertheless, 



Higgs phenomenology may be significantly changed due to mixing of the CP-even and 
CP-odd Higgs states induced by complex mass parameters and coupling constants [|J. 

The Higgs sector of the ^-symmetric NMSSM contains one not removable phase 
which can be chosen to be the phase of Xk*. This phase is not very much constrained 



provided |A| is significantly smaller than one If we allow for Z 3 -symmetry breaking, 
six possibly complex parameters appear in the tree-level Higgs potential. [] Two phases can 
be absorbed by a redefinition of the Higgs and singlet fields. Thus we can take [iB + Ar* 
and kr* being real valued, without loss of generality. We parametrize the Higgs and 
singlet fields according to 

Hl = h x e ie \ H° = h 2 e ie \ S = se i6s , (5.14) 

and define /i = /2e 1 ^, A = Ae"^ A , etc. Furthermore, we introduce 

= 0! +02, = 0i-0 2 . (5.15) 



7 In this section we relax the assumption of universal soft SUSY breaking which would induce corre- 
lations between the various phases at the weak scale. 
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Because of gauge invariance the Higgs potential is independent of the phase combination 
8. Using these definitions the tree- level Higgs potential takes the form 

Vtree = + A V + 2\jjS COs(6 S + <J) X ~ <f>„)) (ft + ft) + \ 2 ftft + 

+2\ks 2 h 1 h 2 cos(# - 28 s + <f) X - <p k ) + 9^-^ft - h 2 ) 2 

o 

+2(Ar* + fiB)hJi 2 cos 8 + 2kr*s 2 cos(20 5 ) + m 2 Hl ft + m 2 H Jil + m 2 s s 2 

2 

+2AA A s/n/i 2 cos(# + 9 S + </>\ + <I>a x ) + ^kA k s s cos(38 s + <f> k + <j>A k ). (5.16) 

Note that in contrast to the MSSM the phase of the //-parameter enters Vtree- Explicit 
CP-violation automatically generates complex vevs for the Higgs and singlet fields. In 
the general NMSSM this is a tree-level effect which should be taken into account. An 
example will be discussed below. EDM constraints on CP-violating phases in the general 
NMSSM depend on the coupling of these phases to the (MS)SM sector. Phases which 
directly enter the squark, chargino and neutralino mass matrices, e.g. M and (8), have to 
obey the MSSM bounds which have been discussed above. The other phases, e.g. 4>k and 
(8s), are less constrained, provided the coupling to the MSSM sector, A, is sufficiently 



small M3 



Even when the Lagrangian is CP-conserving, the CP-symmetry may be spontaneously 
violated by complex scalar vevs, i.e. CP and the SU(2) x U(l) gauge symmetry are broken 
together at the EWPT. In the MSSM radiative corrections to the Higgs potential can 
induce a CP-violating vacuum. However, the required Higgs boson with mass of a few 
GeV is obviously ruled out by experiment |^2[ . The Z 3 -symmetric NMSSM provides only 
limited improvement on the minimal model: radiative corrections generate spontaneous 
CP- violation, provided the lightest neutral Higgs mass is smaller than about 40 GeV |44|] . 
In the general NMSSM the situation is very different. Spontaneous CP-violation (SCPV) 
is possible even at tree-level. However, Higgs spectra consistent with the experimental 
Higgs mass bounds require nearly maximal CP-violation, i.e. (8), (8s) of the order 0(1) 
45|. Phases 8 ~ 0.1 are only possible for Higgs masses smaller than 30 GeV. This finding 



is related to the fact that CP is a discrete symmetry. If CP is spontaneously broken there 
arise two degenerate minima^ with phases ±(#), which are nearby when (8) becomes 
small. From the SM Higgs potential it is known that the squared mass of the Higgs boson 
is proportional to X'(h) 2 , where A' denotes the quartic coupling. In the case of small 
SCPV we can identify \'(h) 2 — > X'v 2 (8) 2 , which up to numerical prefactors of order unity 
is the mass of the emerging light, almost CP-odd Higgs boson. Obviously, for small (8) 
this state becomes arbitrarily light. 

Thus SCPV in the NMSSM does not appear to be very promising. Small phases 
are ruled out because they require a light Higgs boson, while large phases are tightly 
constrained by the EDM experiments. We will see below that this conclusion does not 
apply to SCPV which is only present during the PT. 

We now derive the equations of motion for moduli and phases of the Higgs fields.^ 
Using the definitions (|5.14 ) and ( |5.15| ) we can rewrite the kinetic terms for the Higgs 



8 For that reason explicit CP-violation has also to be present in the model in order to avoid domain 
walls when CP is broken at the EWPT. However, this explicit phase may be much smaller than one. 

9 The singlet field is most conveniently split into real and imaginary part which can be treated along 
the lines discussed in the previous section. 
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bosons in the Lagrangian according to 



C = djixd^fn + dJuVh + ^±^(0^9 + d„9d»9) + f^-J^d^O + .... (5.17) 

Since dgVr = the Euler-Lagrange equation for 9 implies 

(hi + hl)d»e + (hi - h 2 2 )d»9 = c% = const. (5.18) 

In general, the integration constant Cg cannot be set to zero without introducing a pure 
gauge field. However, a non-vanishing Cg would enhance the energy of a bubble configu- 
ration, so we will set Cg = in the following. After elimination of dj) by help of ( 5.18Q 
the equations of motion for hi, h 2 and 9 take the form 



d_ 

dz 



2h 2 h 2 d_ Q 
h\ + h\ dz 



~ ^qVt = (5.21) 



where we restricted ourselves to the case of the static domain wall perpendicular to the 
^-direction ( |4.11|) . The equations for the critical bubble follow along the same lines. In 
eqs. ( |5.19| ) - ( |5.21| ) the phase 9 is dynamical only if both moduli, hi and J12, are different 
from zero. This demonstrates that in the bubble wall 9 cannot be assigned to one of the 
Higgs fields. 

In the previous section we noted that tan/5 varies only slightly in the bubble wall. 
Neglecting this small perturbation the equations ( 5.19Q - ( |5.21|) can be further reduced 



and we obtain for the variation of the phases of the Higgs fields H® and 

59i = sm 2 ((3 T )69, 59 2 = cos 2 ((3 T )59. (5.22) 

Here tan (3t represents the ratio of the Higgs field moduli in the broken minimum at the 
critical temperature. In the limit of large tan(/3r) the variation of 9 is almost completely 
absorbed by 61, while #2 remains more or less constant. For example, in case of tan 0t = 5 
we obtain 59 2 = 0.03850, which is small even for 59 = 0(1). In section 6 we will discuss 
that the variations of complex phases fuel baryon production, rather than the complex 
phases themselves. As a consequence, baryon number generation resulting from a varying 
phase in the top quark mass, m t = y t sin(j3 T )he te2 , will turn out to be highly suppressed 
in the limit of large tan(/3y). Since in the region of the parameter space considered in this 
work the Higgs vev ratio changes only by a few percent when T is raised from zero to the 
critical temperature, this is the case in the regime of large tan /3. 

We now discuss the implications of the general equations of motion derived above. Let 
us start with explicit CP-violation. In the MSSM the complex phases which are induced 
in the Higgs vevs by CP-violating couplings are completely negligible for realistic Higgs 
masses. In ref. || we found this behavior confirmed at finite temperature: 9 < (9(10~ 3 ), 
even for explicit phases of the order of one. The variation of 9 in the bubble is of the 
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Figure 6: (a): Variation of the CP- violating fields (a) and c (b) in the bubble wall as a 
function of h = \/hf + /z| at the critical temperature for M = 0.1, = —100 GeV and 
M = 125 GeV. The remaining parameters are chosen as in fig. 0a. (Units in GeV.) 



same order of magnitude. We conclude that in the MSSM with explicit CP-violation, the 
Higgs vevs can be taken real to very good approximation. 

In the NMSSM the situation is different, because explicit CP-violation is possible 
in the tree-level Higgs potential ( |5.16 ). In general, the behavior of the complex valued 
Higgs and singlet fields in the bubble wall can only be reliably determined by using the 
minimization techniques of ref. |33|]. In order to make this possible, the equations of 
motion which describe s=Re(S'), c=Im(S'), hi, h 2 and 9 have to be included in the T . 
In ref. || this approach was applied to compute CP-violating bubble wall profiles in the 
MSSM and NMSSM. It turned out that also the phase 9 is rather accurately described by 
a kink-ansatz ( |4.12|) . Moreover, it was verified that in the calculation of the wall shape 
small CP-violating phases 9 < O{l0~ l ) can indeed be treated as perturbations of a CP- 
conserving solution. In this case the bubble wall profile can be conveniently determined 
in two steps: first, one computes the CP-conserving profile of the fields hi, h 2 and s, 
where the CP-violating fields 9 and c are set to zero. In a second step, the CP-violating 
constituents of the bubble wall may be computed by simply minimizing the effective 
potential Vr{hi, h 2 , 9, s, c) with respect to 9 and c, while keeping hi, h 2 and s fixed. The 
variation of the Higgs vev ratio can be neglected in this calculation. We emphasize that 
this prescription works only in case of small CP-violation in the Higgs and singlet fields. 

If Vr(hi, h 2 , 9=0, s, c=0) has a smooth ridge, the profile of the CP-conserving fields 
can be obtained from the improved shooting method discussed in section 4. Hence the 
complete CP-violating bubble profile can be computed without minimizing the functional 
T. 

In fig. |] we display an example of a bubble wall in presence of explicit CP-violation 
using <ftn = 0.1, Aq = —100 GeV and Mq = 125 GeV, the remaining parameters are 
chosen as in fig. ^|a. The lightest Higgs boson has a mass = 86 GeV. We find T c = 101 
GeV and v c /T c = 1.77, so the PT is strongly first order. For the wall thickness we obtain 
L w = 5/T. The CP-even singlet field evolves from s = 115 GeV in the symmetric phase 
to s = —79 GeV in the Higgs phase. From fig. ^|a we take 9 = 0.0182 in the broken phase, 
and 59 = 0.0065. In the MSSM, 0^ ~ 1 would be required to induce an effect of the same 
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Figure 7: Transit ionally CP-violating bubble wall profile from the parameter set of table 
[I]. The position variable z is given in units of GeV -1 . 



magnitude. Also the CP-odd component of the singlet field, c, varies in the bubble wall, 
as is shown in fig. [|b. In the symmetric phase we have c = 0, since 0^ couples to the 
singlet only in case of non-vanishing h\^- 

We now turn to the case of spontaneous CP-violation. Although this scenario is 
disfavored at zero temperature, the Higgs and singlet fields may still acquire large complex 
phases at high temperatures, or even only during the EWPT. Since it will turn out in 
section 6 that explicit CP-violation can account for the observed baryon asymmetry only 
in special regions of the NMSSM parameter space, spontaneous CP-violation at finite 
temperature becomes an interesting alternative scenario. 

In the ^-symmetric NMSSM spontaneous CP-violation at finite temperature has 



already been investigated some years ago in ref. [46 , where an effective potential truncated 



at the renormalizable level was used. It was found that spontaneous CP-violation cannot 
occur in the finite temperature broken phase while having a viable temperature zero 
phenomenology, since thermal corrections to the effective potential have the tendency 
to restore symmetries. The authors suggested that spontaneous CP-violation may take 
place during the EWPT ("transitional CP-violation"), i.e. the Higgs and singlet fields 
may acquire complex phases in the bubble wall, although no definite example was given. 
However, the truncated effective potential used in their analysis proved to give misleading 
results in case of the MSSM [|J, as will be discussed below. 

In our analysis of finite temperature spontaneous CP-violation in the general NMSSM 
we use the full 1-loop effective potential ( |3.8| ) without any truncations. It turns out that 
spontaneous CP- violation does not occur for the universal pattern of SUSY breaking dis- 
cussed in section 2, which induces too large masses for the CP-odd Higgs bosons. However, 
spontaneous CP-violation can take place if universality is violated in the singlet sector. In 
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Table 1: Parameter set used in fig. [7] (all dimensionful parameters in GeV). 
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this case A\, and m| are free parameters.^ Fig. [7], which is taken from ref. ||, shows 
an example of a transitionally CP-violating bubble wall for the parameters given in table 
|I|, which was computed using the minimization algorithm of ref. . More precisely, the 
scenario corresponds to a phase transition between an electroweak symmetric high tem- 
perature phase, where CP is violated by a complex singlet vev, and an electroweak bro- 
ken, CP-conserving low temperature phase, i.e. (s, c) sym = (50GeV, 99GeV) — > (s, c)brk = 
( — 122GeV, 0). The critical temperature is found to be 101 GeV, and we have a strong 
first order PT with v c /T c = 1.6. The mass of the lightest CP-even Higgs boson is 84 GeV 
which because of its reduced coupling to the Z-boson is still compatible with the experi- 
mental data. There is also a rather light CP-odd state with a mass of 105 GeV, which is 
almost a pure singlet. The bubble wall is found to be rather thin with L w ~ 3/T. In our 
example spontaneous CP-violation is triggered in the singlet sector. The small coupling 
A communicates it to the Higgs fields. As a consequence, the almost maximal phase of 
the singlet field, 9 S , induces only a small phase 9 ~ 1/20 ~ \9 S in the Higgs sector. Up to 
now we only identified one particular region in the NMSSM parameter space where transi- 
tional CP-violation occurs. A more systematic study would be very desirable. Especially 
it would be interesting to find an example where transitional CP-violation originates from 
the Higgs fields rather than the singlet field. In that case 9 is expected to be of the order 
one. 

Let us finally comment on transitional CP-violation in the MSSM, where there also 
have been claims in the literature that transitional may occur j|7|]. In ref. |48| even an 
example of a transitionally CP-violating parameter set has been presented. All these 
studies were carried out using an effective potential that has been truncated at the renor- 
malizable level, i.e. the logarithmic terms, etc., were neglected. In ref. |5[] we performed a 
systematic search for transitional CP-violation in the MSSM using the full 1-loop thermal 
effective potential without finding any viable parameter set. In particular, the example 
of transitional CP-violation presented in ref. could not be confirmed. This discrep- 
ancy appears to be due to the different approximations which have been made concerning 
the effective potential. Our results indicate that the truncated effective potential is not 
appropriate to discuss transitional CP-violation. 



6 Baryon Asymmetry in the Semi-Classical Limit 

6.1 Generalities 

Between the initial nucleation and the completion of the phase transition expanding bub- 
bles convert the symmetric phase into the broken phase. The early universe is filled with a 
hot plasma of particles, most of them have rather different masses and mixings inside and 
outside the bubble. Thus the bubble wall behaves like a potential on which the particles 
scatter. At the bubble surfaces the plasma is thrown out of equilibrium by the motion of 
the phase boundary. The higher pressure inside the bubble tends to accelerate the wall, 
while the interaction between the wall and the particles in the plasma dissipates energy 
and slows the wall down. Finally, a stationary situation is reached, where the different 
forces balance and the wall propagates with constant velocity v w . In general, the wall 
velocity depends on the shape of the effective potential and on the composition of the 

10 It is even sufficient to violate universality only via m|. 
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plasma. 

In the symmetric phase baryon number violation occurs frequently due to hot sphaleron 
processes. A baryon asymmetry is produced, if the departure from equilibrium at the 
phase boundary biases the baryon number violating processes in a CP-violating fashion. 
The most efficient mechanisms of baryon number generation rely on transport. If CP is 
violated in the interaction between the bubble wall and the particles in the plasma, dif- 
ferent population densities for particles and anti-particles are induced in front of the wall. 
The difference between the particle and anti-particle populations is then transported into 
the symmetric phase, where it biases rapid baryon number violation. The total amount 
of baryon number which gets produced during the phase transition depends crucially on 
the shape and motion of the bubble wall. 

Different methods have been suggested in the literature in order to describe the effects 
of CP-violating interactions between particles in the plasma and the propagating bubble 
wall, which finally generate CP-violating source terms that enter the diffusion equations 
for particle transport. Depending on the properties of the bubble wall, in the first place 
its velocity v w and thickness L w , different approximations may be applied. The most 
rigorous approach is based on the closed time-path (CTP) formulation of non-equilibrium 
quantum field theory [|5C]. It leads to a set of quantum Boltzmann equations [|51] which 
describe the temporal evolution of particle densities including particle number changing 
interactions and CP- violating source terms |52 |. However, the approximations and results 
are still controversely discussed. There have also been efforts to use the CTP formalism 
only in calculating the CP-violating source terms which in turn are inserted into classical 



Boltzmann equations [49 



If the thickness of the bubble wall L w is smaller than the mean free path I of the 
particle under consideration ("thin wall"), one may neglect the influence of the plasma 
during the scattering of the particle off the wall. The interaction of a fermion or boson with 
the wall can then approximately be described by using a "free" Dirac or Klein-Gordon 
equation, respectively. CP-violation is encoded in different reflection and transmission 
coefficients for particles and anti-particles [|53|-|57j . In ref . [|5£|] also the effects of thermal 



scattering were taken into account, which due to decoherence have a negative impact on 



the generation of a CP- violating observable fo9fl . Baryon production is fueled by a net 
flux of charge into the symmetric phase ("charge transport mechanism" |[53||), which is 
induced by the reflection asymmetry. This current is subsequently inserted into a set of 
classical Boltzmann equations which determines the evolution of the particle distributions. 
Recently, reflection and transmission probabilities have also been calculated by using 
standard quantum field theory methods |p0[ . There, thermal scattering, i.e. damping, 
was included via an imaginary part of the self-energy. 

As the bubble wall becomes thicker, L w ~ I, interactions with the plasma must in- 
evitably be taken into account, when a particle encounters the propagating wall. In 
general, this requires non-equilibrium quantum field theory methods [f^ . However, if 
L w ^> 1/T, most particles have inverse momenta 1/p <C L w and may therefore be treated 
semi-classically pT|-|5^]. Using the WKB approximation one can derive dispersion rela- 
tions which in the presence of CP-violation are different for particles and anti-particles. 
The dispersion relations then enter classical Boltzmann equations, leading to a unified 
description of CP-violating source terms, particle scattering and transport. Recently, it 
was clarified how the semi-classical description arises from the quantum Boltzmann equa- 
tions in the limit of a slowly varying background field (thick wall regime) ||65||. Of course, 
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even in the thin wall limit the high-momentum particles may be described semi-classically. 
However, in that case they only give a sub-leading contribution to the CP-violating source, 
which in this case is dominated by the low-momentum particles. Even in case of thick 
walls there may be important contributions from low-momentum modes neglected in the 



semi-classical treatment, which according to ref. [[52| significantly enhance the produced 
amount of baryon number. There is also a very recent calculation [66| showing that gauge 
fields in the hot phase lower significantly the wall velocity, such also leading to a more 
efficient baryon production. In ref. |67j the semi-classical approximation to particle dy- 



namics has been applied to calculate the bubble wall velocity in the SM. In ref. |68| it 
was shown that the stops in the MSSM considerably lower the wall velocity. 

Notice that different particle species with different (gauge) interaction have different 
mean free paths. Thus whether one is in the thin or thick wall regime depends on the 
particle species under consideration. In the (N)MSSM especially particles with strong 
interactions (quarks and squarks) require a thick wall treatment. 

There is some disagreement among the groups that have estimated the baryon asym- 



metry generated during the EWPT in the MSSM. In refs. j|9], [5^, the CP- violating 
source terms are proportional to the variation of tan/? in the wall, which according to 
the discussion of section 4 causes a suppression of at least C(10~ 2 ). This dependence 
is due to taking into account only the leading order of the Higgs insertion expansion, 
which has been used to calculate source terms. At higher orders in the expansion there 
arise contributions, which escape the 5/3 suppression, as was recently shown in ref. ||60|| . 
However, these corrections turn out very small and are competitive to the leading order 



contributions only in case of 5(3 < O(10 ). In the recent paper [71] it was argued con- 
vincingly that the (tan 0)' source cancels and that one has to consider closely a source 
term symmetrical in the two Higgs fields. 



In this section we generalize the method introduced in ref. |6J] to calculate the baryon 
asymmetry in the NMSSM, which has not been considered so far. Different from ref. [|64[ 
our formulas also cover the case of CP-violating bubble walls. We obtain a non-vanishing 
baryon asymmetry even in case of constant tan f3, which is due to the variation of the 
singlet field or the presence of CP-violating bubble walls. 

6.2 WKB approximation and dispersion relations 

The idea of deriving CP-violating source terms using a semi-classical approximation was 
developed in context of the two Higgs doublet model |k| and afterwards applied to the 



charginos in the MSSM In the following we review the derivation of semi-classical 
dispersion relations and generalize this method to account for the effects of CP-violating 
bubble walls, which potentially are present in the NMSSM. These additional sources of 
CP-violation will turn out to be very helpful in order to generate a baryon to entropy 
ratio in the observed range: r\ B = tib/s ~ 2 — 7 x 1CT 11 p9| . We also give the dispersion 
relations in the bosonic case. In the next section these dispersion relations enter the 
classical Boltzmann equations that describe particle transport induced by the propagating 
bubble wall. 
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6.2.1 The fermionic case 



We start the derivation of fermionic dispersion relations in presence of CP- violation with 
a simple example: A single (Dirac-) fermion = ij] ai x a ) that couples only to one of 
the two Higgs doublets, which we denote by H. The most prominent realizations are the 
quarks and leptons of the (N)MSSM in the presence of a CP-violating bubble wall. Due 
to its coupling to the Higgs, y, the fermion obtains a mass proportional to the Higgs vev, 
M = yH. During the passage of the bubble wall the fermion mass becomes space-time 
dependent. Assuming that the bubble has grown to macroscopic size and reached its 
final velocity, we can neglect the curvature of the wall and boost to the rest frame of the 
bubble. Then the fermion mass only depends on one position coordinate, which we denote 
by z, i.e. M = M(z). As a consequence, the energy of the particle, E, and its momentum 
perpendicular to the ^-direction, p±, are constants of motion. Writing M = me , possible 
CP-violation is encoded in a non-vanishing phase Any constant 6 can be absorbed 
by a redefinition of the fermion field. When couplings to additional fields are neglected 
only a varying phase contributes to CP-violation. 

In the presence of a complex mass term the fermion field is described by the free Dirac 
equation 

(if^ - P L M - P R M*)* D ee ( ia _M* ) (?) = °- ( 6 - 23 ) 

Notice that at this level all interactions between the particle and the plasma are neglected. 
The scattering effects will be accounted for in the next section when Boltzmann equations 
are written down that describe the local phase space distributions. We are working in the 
chiral representation of the 7-matrices. Exploiting conservation of energy and boosting to 
the Lorentz frame where p± = we can take the ansatz ^ d = e~ tEt ^{z) and are left with 
a one dimensional problem. The interaction between the fermion and the wall conserves 
the z-component of the spin, S z . Thus eq. ( p.23j ) splits into two equations 



id& = ±Q(z)t ± , Q(z) = [ ) ( 6 .24) 



where £+ = (£1, £3) and = (£2, £4) are the S z = ±| components of £. To solve eq. ( |6.24|) 



one brings the ^-dependent matrix Q{z) = D(z)Qe>(z)D(z) 1 into a diagonal form, where 

HI 



, VE 2 - m 2 \ / coshX e-' e sinhX , 

Ql) " [ -VE 2 -m 2 h { e ie smhX cosh A ' {b ' 2b) 



and tanh2X = m/E. In the local helicity basis, £± = D 1 £-j_, the Dirac equation ( |6.24| ) 
takes the form 

ihdj ± = (±Q D - D-Hhd z D)l ± (6.26) 

which still is an exact equation. In general, the correction term D~ 1 ihd z D caused by 
the position dependent field redefinition is not of diagonal form. The two components 



11 Here 9 denotes the phase of the Higgs boson the fermion is coupling to. It should not be mixed up 
with the common phase of both Higgs fields defined in eq. ( |5.15 ) which is denoted by the same symbol. 
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of £± are still coupled. However, the off-diagonal part is proportional to d z D ~ D / L w . 
Typical momenta of the particles in the plasma are of the order of the temperature T, 
which is much larger than 1 / L w for the bubbles under consideration. We therefore expand 
eq. ( |6.26| ) in powers of d z or more precisely in powers of h (WKB approximation) that we 
already reintroduced for that reason. 

To order (h)° we can neglect the D~ 1 hid z D contribution. Thus the two components 
of £± decouple in eq. ( 6.26|) . Inserting the WKB ansatz for the fermion field 



|£) = ( iV-tS'toVW, g 2 ) = (°\ e -U'p.{*W ( 6 .27) 



S>J U 

into (|6.26| ), we obtain the dispersion relations p z (E) 

i?,^: Vz = VE 2 - m 2 



C+M-^ Vz = -VE 2 - m 2 . (6.28) 



The momenta p z are the eigenvalues of the matrix entering the RHS of eq. ( J6.26 ). The 



eigenfunctions may more transparently be labeled by the chirality states they correspond 
to in the limit m — > 0, e.g. ~ £i ~ etc. The dispersion relations may then be com- 
bined to p z = sgn(p 2 )V-E' 2 — which holds for left-handed particles 77 and right-handed 
particles % Obviously, in the classical limit the ^-dependence completely disappears, 
demonstrating that CP-violation is indeed a quantum-mechanical phenomenon. 

To solve the Dirac equation (|6.26|) to order h we have to take into account the 
D~ 1 hid z D term which reintroduces a coupling between the two components of £±. The dis- 
persion relations p z (E) are obtained from the eigenvalues of the matrix ±Q d — D~ x ihd z D . 
Since to order h the off-diagonal terms do not contribute, we are left with 



L (77) : p z = sgn(p z )V E 2 — m 2 — ^'sinh 2 X, 

R (x) : p z = sgn(p z )VE 2 - m 2 + W sinh 2 X, (6.29) 



where 8' = d 7 9 and 



, , , E — VE 2 - m 2 



sinh 2 X = / (6.30) 

2V^ 2 - m 2 V ; 

Again the states are labeled by their asymptotic chirality properties. Notice that the vari- 
ation of m, which is encoded in d z X, drops in the dispersion relations. The CP-violating 
part of the dispersion relation, Ap z = W sinh 2 X, is proportional to the derivative of the 
phase 6. Thus only a varying phase contributes to CP-violation in the semi-classical limit. 
Furthermore, CP-violation is proportional to sinh X, which guarantees that its effect is 
turned off in the limit m — > 0, where 9 is no longer well defined. Because of the different 
dispersion relations, left- and right-handed particles feel a different (semi-classical) force 
in their interaction with the wall. For the anti-particles, fj (L) and x (R), the CP- violating 
part comes with the opposite sign. Besides the force term there is a second manifestation 
of CP-violation: The phase 9 enters also the transformation matrix to the helicity basis 
D. Since the interaction eigenstates are different from the helicity states, particles and 
anti-particles interact in a different way with the surrounding plasma. This generates a 



CP- violating source term which drives "spontaneous" baryogenesis [70 
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As pointed out recently in ref. [71], one should better use the kinetic momentum 
Pkin = ^'Wgroup = m ff instead of the canonical momentum p (which we used up to now) 
in the quasi-classical limit of particles in the Boltzmann transport equations. This is quite 
in the spirit of the correspondence principle of basic quantum mechanics. Calculating the 
(inverse) group velocity from ( |6.29|) , and using (|6.30|) , the kinetic moment beyond the 
zeroth order contains (order K) correction terms 

Apkin = ± ^ . (6.31) 

In the following we will need the dispersion relation for energy in terms of momentum 
which to order h is obtained as 

9 

m 



E ± = E ±AE= ^ in + m 2 ± signal? +m2) (6.32) 

E + is the energy of left-handed particles and anti-particles, whereas E_ corresponds to 
right-handed particles and anti-particles.Q In the derivation of eq. ( |6.32| ) we transformed 
to a general Lorentz frame with non-zero momentum parallel to the wall. When boosting 
to the plasma frame the dispersion relation ( |6.32|) is preserved to linear order in v w and 

In our derivation of the dispersion relations we followed basically ref. . CP- violation 
arises because of a position dependent phase in the transformation from the interaction 
(i.e. chirality) states to the local mass (i.e. helicity) eigenstates. In ref. f63| a slightly 
different approach was used to obtain the dispersion relation of a fermion (e.g. top quark) 
in the presence of a CP-violating Higgs field background, H, in context of the 2HD 
model: The complex phase in the fermion mass was removed by a gauge transformation, 
which induced a gauge field in the kinetic term of the fermion. The dispersion relation 
obtained with this technique is completely analogous to our result ( ft.29 ). However, it is 



not clear how to generalize this method to cover the case where several species mix, as 
occurs with the gauginos and Higgsinos in the (N)MSSM, or where the NMSSM singlet 
field background contributes to CP-violation. In both situations the CP-violating phase 
cannot be removed by a gauge transformation. On the other hand, the method described 
above is still applicable [B3|. 



We are now in the position to address the problem of mixing Dirac fermions *&di, 
where the "flavor" index I — 1, N. In order to solve the corresponding Dirac equation 
we diagonalize the matrix 

OM-U, -£ W ) <«•»> 

which is a straight forward generalization of (|6.24|) . Here 1 denotes the unity matrix in 
flavor space and M is a position dependent complex mass matrix for the fermions, without 
further specifications. As a first step we write M and in terms of two unitary matrices 
U and V, and a diagonal matrix = me 10 

M = VMoU 1 , M" 1 " = UM^Vt (6.34) 



12 



Notice that in ref. [p4| the phase 9 has been defined with the opposite sign. 
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U and V may be obtained by diagonalizing the hermitian matrices M^M and MM', 
respectively. The transformation to the helicity basis takes the form 



coshX -e- w sinh X \ / TJt 
-e i0 sinhX coshX j { V f 



T-^ ± (6.35) 



where X is a diagonal matrix in flavor space, obeying tanh(2X) = m/E. Like in the single 
fermion case the CP-violating part of the dispersion relation is encoded in the diagonal 
elements of T~ l ihd z T. We find 



i { ± } : Pz = ±^E 2 -mj- 9\ sinh 2 X l - cosh 2 X I [\jHd z \J] I + sinh 2 Xi[VUd z V]i, 

lf> : Pz = T JE 2 -mj + 9j sinh 2 X/ - cosh 2 X^V^V], + sinh 2 X / [U t i9,U] / {6.36) 



Here the subscript I on the RHS denotes the 1th diagonal element of the corresponding 
matrix, i.e. mi = m/j, and we dropped the factor h. In addition to the 9' contribution, 
which we already encountered in the case of a single Dirac fermion, the dispersion relations 
receive corrections due to the position dependent rotations in flavor space, U and V. 
In eqs. ( p. 35 ) and (|6.36|) we allowed for complex phases in the diagonal matrix Md. 



Since U and V contain 2X 2 — 2 real parameters it would be sufficient to take real values 
in ~M.£> in order to reproduce the 2X 2 real parameters of M. However, we have to omit 
the N — 1 rotations in U and V which belong to the Abelian subgroup of SU(N). These 
are not related to CP-violation, but rather correspond to artificial redefinitions of the 
interaction (chirality) states, as can be deduced from the case vanishing fermion mixing. 
Thus allowing for complex values for the masses in we end up with a one by one 
correspondence between the X 2 parameters in M and the ones contained in U, V and 
M D . 

Having resolved the ambiguity in the definition of the helicity basis, it is possible to 
determine the matrices U, V and Mo numerically. The dispersion relations are then 
easily obtained by evaluating the expressions ( |6.36[ j . Here we will concentrate on the 
charged Winos and Higgsinos in the (N)MSSM. Since they mix via the 2x2 chargino 
mass matrix 

C =- + <? W > h ^{g 2 (H°y ^ + XS ){hi) (6 ' 37) 

analytic formulas can be obtained. The charged Winos and Higgsinos can be combined 
to a Dirac spinor ^ D = (^ L ,^ R ) T , where ^ L = (W + , h^) T and ^ R = (W~,iq) T . We 
parametrize the SU(2) matrices U and V by 

TT / cos a e~ n sina \ _ r / cos b e~ l& smb \ ,.„„> 
U = j~ . , V = iS . (6.38) 

V -e'sma cos a J \ —e smb coso I 

According to the discussion above we dismissed phases multiplying the cos-terms. The 
diagonal elements of \J^id z U which enter the dispersion relations (|6.36|) read 

[U f ^ 2 U]i = -[\] ] id z \J] 2 = -7 sin 2 a. (6.39) 

We observe that only the derivative of the complex phase 7 contributes, whereas the 
derivative of a drops. Similar relations hold for "VHdzV. Inserting these expressions 
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into eq. (|6.36|) we obtain the dispersion relations for left- and right-handed particles and 
anti-particles 



p z = sgn(p z ) J E 2 — m 2 — (6'j + 8' sin b) sinh Xj + 7'sin a cosh Xj, 



Lj : p z = sgn(p 2 ) \ E 2 — mj + (9'j + 8' sin 2 b) sinh 2 Xj — 7' sin 2 a cosh 2 Xj, 



E 2 


2 

— mj 


E 2 


2 

— mj 


E 2 


2 

— mj 


'e 2 


2 

— mj 



In the symmetric phase L 2 and R 2 evolve to the left-handed Higgsinos states h 2 and h^, 
respectively. The flavor transformations U and V are related to the parameters of the 
chargino mass matrix (|6.37| ): 

sin 2 a = 2|v4| 2 /A(A + A) with 

A = g 2 ((M 2 H 2 )* + (n + \S)H° 1 ) 

A = \M 2 \ 2 -\^ + XS\ 2 + g 2 2 (\H 1 \ 2 -\H° 2 \") 

A = (A 2 + 4|A| 2 ) 1/2 (6.41) 



and 7 = argA This gives 

i sin 2 a = 2lm(A*A') /A(A + A) 

and there are similar relations for sin 2 b, 5 and 8' sin 2 b exchanging a and b, H® and H® , 7 
and —8. The mass eigenvalues read (in non-symmetric notation) 

r,^ i ,^cosa / TT nN*sina 

[M D ] n = M 2 --g 2 (H°y -e~ l \ 

cos b cos 

[M D } 22 = ^ + \S) C ^ + g 2 (H?y S ^-e^. (6.42) 

coso coso 



For A = these expressions agree with the results of ref . |7T[ ■ In this case one is left with 
the MSSM, and the phases 7 and 8 only vary due to a change in the Higgs vev ratio tan (3 
or because of transitional CP- violation in the bubble wall. The first contribution is highly 
suppressed, since the variation of (3 is at most ~ 1CT 2 for realistic Higgs masses p3|,|35||, 
while transitional CP- violation most probably does not occur at all in the MSSM ||. On 
the other hand, the contribution to the chargino dispersion relations stemming from the 
variation of the complex phases in requires only explicit CP-violating phases in \i 
or M 2 . Eq. ( |6.42j ) demonstrates that even though the phases in the two terms entering 
[Md] n t 22 are position independent, their contribution to the resulting phase varies due 



to the change in the (real) Higgs vevs 64]. According to the discussion in section 5, the 
corrections due to complex Higgs vevs induced by the explicitly CP-violating phases can 
also be safely neglected in the MSSM p. 

In the NMSSM CP-violation enters the dispersion relations in several new ways. As 
already discussed in section 5, spontaneous CP-violation in the bubble wall occurs for 
specific values of the SUSY parameters, which leads to a variation of 7, 8 and 6\ y2 . If this 
effect is present, it dominates the CP- violating part of the dispersion relations. Even in 
the absence of transitional CP- violation and with constant tan (3 there are contributions 
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Figure 8: CP- violating contributions to the chargino dispersion relation for the explicitly 
CP- violating parameter set of fig. || (a) as a function of z for p = T c and (b) as a function 
of p for z = 0. The dashed-dotted, dashed and solid lines represent the helicity and flavor 
contributions ( |6.40|) and the result using the kinetic momentum (|6.43| ), respectively. 



to 7' and 5': The complex phases induced in the vevs by the explicitly CP- violating phases 
may considerably change in the bubble wall according to section 5. Also the phase of the 
effective /i-term, ju + XS, is position dependent because of the variation of the singlet vev. 
However, this effect is suppressed when the coupling A is small. 

Again we calculate the (inverse) group velocity, now from ( |6.40| ). ^-independent 



terms drop out, i.e. cosh 2 Xj can be substituted by sinh 2 X/. The kinetic momenta 
beyond the zeroth order ( |6.29D then contain (order h) correction terms ±(0/ + 5' sin 2 b — 



7' sin 2 a)m 2 /2E(E 2 — m 2 ) 1 / 2 . The dispersion relation can be inverted, and to leading order 
in the derivatives the CP-violating part of the dispersion relation for the eigenstate L 2 
corresponding to h\ in the symmetric phase is 



AE = -sign( p z )(6' 2 + 5' sin 2 6-7 sin 2 a) , ™ 2 — ^ (6.43) 

2 (Pkin + m 2) 



where m 2 = IfiWzj^l 2 from eq. (|6.42|) . For pkin rn this is 2\p z \/p times the result one 
would obtain with canonical momentum after the substitution of cosh 2 by sinh 2 in ( |6.40| ) 
(see fig. |8|b). Note that AE is now totally symmetric under the exchange of Hi and H 2 . 
This destroys the most prominent source term ~ HiH' 2 — H[H 2 of older work. 

In fig. § we compare the CP- violating contributions to the chargino dispersion relations 
stemming from the local helicity and flavor transformations ( |6.40| ) with the result from 
using kinetic variables ( |6.43|) . We assume that p± = 0, and approximate the bubble 
wall profile by a kink-ansatz with a common wall thickness for all fields, i.e. we take the 
straight connection between the symmetric and the broken minimum in field space. In 
this approximation tan (3 is automatically constant along the wall. We study the case of 
explicit CP- violation induced by ^ = 0.1, which has already been considered in fig. |6] in 
section 5. It is characterized by L w = 5/T c and T c = 101 GeV. Fig. ||a shows the position 
dependence of AE for p = p z = T c . In fig. |8|b we present AE as a function of momentum 
for z — 0. One observes that AE is dominated by the helicity contribution proportional 
to 8' 2 (dashed-dotted line), while the flavor contribution (dashed line) only gives a small 
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correction. The result from the kinetic approach is larger by a factor between 1 and 2. 
For the transitionally CP-violating parameter set of fig. [7| the results are similar. The 
amplitude of AE, however, is by a factor of about 7 larger, due to the larger amount of 
CP- violation and the thinner bubble wall (L w = 3/T). In any case, AE provides only a 
small correction to dispersion relation. 

We close this paragraph by briefly considering the case of Majorana fermions. Ma- 
jorana spinors have the special property of being invariant under charge conjugation, 
i.e. particles and anti-particles are contained in the same four-component spinor. The 
mass matrix of Majorana particles is symmetric, and the entries are complex numbers, in 
general. In the NMSSM, the neutralinos are Majorana particles with a 5 x 5 mass matrix. 
From the symmetry property of the mass matrix we conclude MM' = MM* = (M*M)*. 
The flavor transformations therefore obey V = U*. Inserting this result into the gen- 
eral dispersions relations ( |6.36j ) (or into its two-dimensional version ( p.40| )) we find that 
"left-handed" and "right-handed" states acquire exactly the opposite CP-violating con- 
tribution, in contrast to the general Dirac case. This is not surprising, since in case of 
Majorana fermions the two helicity states describe particles and anti-particles, respec- 
tively. 

6.2.2 The bosonic case 

In supersymmetric models the scalar superpartners of the top quarks may give impor- 
tant contributions to the CP-violating source that fuels baryon production. The stops 
efficiently interact with the bubble wall via the large top Yukawa coupling. They contain 
many degrees of freedom, and due to renormalization group flow the right-handed stop 
is probably lighter than any other squark. In the following we consider the interaction 
between the squarks and the bubble wall in the semi-classical limit. Our approximation 
is valid if the thickness of the bubble wall is much larger than the typical wavelength 
of particles in the plasma which is of order 1/T. As in the fermionic case, we derive 
dispersion relations in which the CP-violating part is generated by varying phases in the 
scalar mass matrix. 

Consider N complex scalar fields A = (A%, ...A^) T with mass term A f M 2 A. The 
entries of the hermitian mass matrix M 2 may be complex. Due to the interaction with 
the wall, M 2 depends on the position variable z. We represent M 2 by a unitary matrix 
U and a real diagonal matrix M^ = diag(M 2 ?1 , Mp N ) 

M 2 = UM^U f . (6.44) 

Like in the fermionic case we omit the N — 1 rotations in U which belong to the Abelian 
subgroup of SU(N). These transformations artificially redefine the interaction eigenstates. 
Their phases are not related to CP- violation, as can be deduced from the case of vanishing 
mixing. We then have a one to one correspondence between the N 2 real parameters in 
M 2 , and in U and M^. 

From the Klein-Gordon equation for the scalar fields (f^d^d^ + M 2 )A = we obtain 
the equation of motion for the local mass eigenstates A = U'A 

[ - h 2 {d 2 z + 2\J ] d z \Jd z + Ut^U) + M 2 - E 2 ] A(z) = 0. (6.45) 

To derive this equation we implemented energy conservation in the wall frame by the 
ansatz A(t, z) = A.(z)e~ lEt and boosted to the Lorentz frame where there is no momentum 
parallel to the wall, i.e. p± = 0. 
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We solve eq. ( |6.45|) by the WKB method. To order O(ftP) we can neglect the Ut<9 2 U 
and XJ^d z XJ contributions and find the wave functions A^(z) = e/ exp(— i f z p z (z')dz'), 
where is the 7th unit vector in flavor space. Furthermore, we obtain the dispersion 
relation p 2 = E 2 — M^j. As we already found in the fermionic case, CP- violation vanishes 
in the classical limit. 

We now include the order 0{h 1 ) corrections. We still can neglect the U^c^U term, 
while the XJ^d z XJ contribution becomes relevant now. Its off-diagonal entries reintroduce 
a coupling between the N components of eq. (|6.45|) . In the dispersion relation only the 
diagonal part of XJ^dzXJ enters. We obtain 

A« : E 2 = M 2 DI +p 2 z + 2p z [UHd z U]i, (6.46) 

where [UHd z U]i denotes the Jth diagonal entry of the corresponding matrix. To transform 
the dispersion relations to a general Lorentz frame with non-vanishing p± one has to 
replace E 2 by E 2 —p\- 

Let us evaluate the general expression (|6.46|) for the case N = 2. Using the represen- 
tation of U from eq. fl6.38| ), we obtain 

A (1 ' 2) : E 2 = M 2 Di 2 +p 2 z - 2p zl ' sin 2 a, (6.47) 

Again, CP-violation arises due to the varying phase 7 in the transformation to the local 
mass eigenstates. A variation in a gives no contribution to the dispersion relations. For 
the anti-particles, A^ 1,2 -', the CP- violating part in eq. ( |6.47| ) enters with the opposite sign. 
If we apply these expressions to the top squarks of the (N)MSSM with mass matrix 

r - (f* t c \ ( m ^ m ^ \( * 
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y t (fi + XS)H 1+ y t A;(H 2 r, (6.48) 



the parameters of U are given similarly to ( 6.41|) by sin 2 a = 2|y4| 2 /A(A + A), 7 = arg A 
now with A = y t ((fi + \S)H° + A* t H$*). Solving (|6~4p for p, one obtains p(E) = (E 2 - 



M^,^) 1 / 2 + 7' sin 2 a to first order in 7'. The group velocity (||) 1 is independent of 7', 
thus stops do not contribute to CP-violation in this order in the kinetic approach. 



6.3 Diffusion equations 
6.3.1 The fluid approximation 

In this section we study the coupled differential equations that describe particle inter- 
actions and transport during the phase transition. We treat the plasma as consisting 
of quasi-classical particles with definite canonical position and momentum. This is an 
approximation to quantum Boltzmann equations which have to be discussed in princi- 
ple. This picture is justified for thick walls (p ^> 1/L W ) if it predicts a sizable effect, 
not dominated by non- leading terms in the derivative expansion [p5| . The dynamics of 
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the particles is then governed by the dispersion relations E(x,p) derived in the previous 
section. 

The information about the particle distributions is encoded in the phase space densities 
fi(x,p,t). Their temporal evolution follows from the Boltzmann equation 

d t fi = (dt + £-d 3 + p- d p )fi = Ci[f}. (6.49) 

The time derivatives of position and momentum obey the Hamilton equations x = 
dfiE(x,p) and p = —dgE{x,p). The Boltzmann equation can in principle be solved nu- 
merically. However, to make it analytically tractable we use the fluid-type truncation 

m 

fi(x,p,t) = em _J z _^ )±l (6.50) 

for the phase space densities of fermions (+) and bosons (-) in the rest frame of the plasma. 
Here Vi and //j denote the velocity perturbations and chemical potentials for each fluid, 
respectively. We also split E^ into a dominant part E 0i = ^Jp 2 + mf and a perturbation 
AEi ~ d z 9 which is related to CP-violation. The chemical potentials are the central 
quantities that finally will determine the baryon asymmetry. The velocity perturbation 
on the other hand, is only introduced to allow the particles to move in response to the 
force, giving rise to chemical potential perturbations. The fluid truncation is valid as long 
as perturbations beyond the ansatz (|6.50|) are attenuated faster than chemical potential 
perturbations. As discussed in ref. |k| this requires v w < where the diffusion constant 
D will be introduced below. 

We are looking for a "stationary" solution of the Boltzmann equation, because at 
late times the wall moves with constant velocity v w . This means that any explicit time 
dependence enters in the combination z = z — v w t. Inserting the fluid ansatz into the 
Boltzmann equation ( |6.49 ) we obtain to linear order in the perturbations AEi, Hi and Vi 

fj - v w (AE> - v' tPz - /<■) - j±Mp M + /4) + ^p^j = d\f\ (6.51) 

where f'± = df±/dEo = —(3e^ E °/(e^ E ° ± l) 2 denotes the derivative of the unperturbed 
Fermi-Dirac or Bose-Einstein distribution. The remaining primes in eq. ( |6.51|) denote d s - 
More precisely, eq. ( |6.51| ) is the difference of the corresponding equations for particles 



and anti-particles. The parameters AE i: , /ij and Vi therefore represent the differences of 
these quantities for particles and anti-particles. For that reason the term f±v w E' 0i , which 
provides the main contribution to friction in the calculation of the wall velocity ||67|| , 



cancels out. In order to obtain differential equations for the perturbations we average 
over momentum, weighting eq. ( |6.51| ) by 1 and p z , respectively, and receive 

2 

-(^-X + k^ = (Q, (6.52) 
v w {plH-(^)^-v w (p z AE' 0i ) = ( Pz Q. (6.53) 

&0i 

The average is defined according to 

J d 3 pf' + {m = 0) J d 3 pf ± 
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The statistical factor k is 1 for massless fermions, 2 for massless bosons and exponentially 
small for particles much heavier than T. In the derivation of eqs. (|6.52 ) and ( |6.53| ) one 



uses that the energy perturbations AEi are odd in p z , according to the results of the 
previous section. Furthermore, we neglected the contribution of the last term in the curly 

2 

brackets of eq. ( |6.51| ) which turns out to be numerically less important than the (J^r)^ 



term taken into account in eq. ( |5.52|) . 

To linear order in the perturbations the collision terms on the RHS of the Boltzmann 
equation take the form 



< c *> = E r "E^> (pA) = v,piJ2 Te P - ( 6 - 55 ) 

p j p 

where T p denotes the rate of the process p. In the sum over chemical potentials, incoming 
particles enter with positive sign, outgoing particles with negative sign. We distinguish 
between elastic interactions with rates T e and decays (inelastic interactions) with rates T d . 
Since elastic scattering conserves the number of particles these processes do not contribute 
to (Ci). In the hot electroweak plasma the elastic processes dominated by gauge boson 
exchange reactions are much more efficient than the particle decays coming from Yukawa 
interactions, sphalerons, etc. It is therefore justified to neglect the contribution of inelastic 
processes to (p z Ci)- 

In the collision integral a delta function 5(^2 Pi) usually represents the conservation 
of energy and momentum in the interactions. However, in the plasma frame this is no 
longer true, because of the moving bubble wall. The Boltzmann equation is formulated 
in terms of mass eigenstates which change in space and time due to the varying Higgs 
(and singlet) vevs. As a result, additional CP-violating contributions to the collision term 
arise, which are related to what has been dubbed "spontaneous" baryogenesis [7D|. The 
final form of the collision term is therefore given by [|64| 

(d) 



E r "(E^ + AE wV (6-56) 
p V j / 



where AE sPtP denotes the bubble wall induced deviation from energy conservation in the 
corresponding process. 

Finally, we reduce the two coupled transport equations fl6.52| ) and ( |6.53| ) to a single one 
by differentiating (|6.53| ) and eliminating Vi in favor of //j. We approximate the thermal 



averages according to (p 2 z /E ) ~ Ki(p 2 z / E ) , (p 2 z ) ~ Ki(p 2 z ) , where the subscript "0' : 
denotes averaging with the massless, unperturbed Fermi-Dirac distribution. Defining the 
diffusion constant D { = Hi(p 2 z / E )l/ \p 2 Tf) we find the diffusion equation [Q 

- Ki(Diix'( + Vvii'i) + r p Yl & = S " 

p j 

Si = j^yiPzAE')' -J2^ AE ^,p)- (6-57) 
\Vzl -fro/o „ 



In order to obtain eqs. Q6.57 ) we neglected derivatives of the CP-conserving thermal 



averages and rates, and left aside ratios of inelastic to elastic scatterings. We also ignored 
terms beyond leading order in the wall velocity, an approximation that certainly breaks 
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down if the wall velocity approaches the speed of sound in the plasma, v s = 1/ ~ 0.58. 
If the wall moves faster than v s , perturbations cannot propagate in the region in front of 
the wall any more. Notice that the first contribution to the CP-violating source term Si 
in ( |6.57| ), which is due to the semi-classical force, is proportional to the diffusion constant. 
This is simply because particles must move in order to build up perturbations. The second, 
"spontaneous" source term is independent of transport properties of the corresponding 
particles and therefore dominates in the limit of inefficient transport. In Si the thermal 
averages over the CP- violating energy perturbations AE are performed using the massive 
distribution functions in order to account for Boltzmann suppression of heavy particles. 



6.3.2 Diffusion equations for supersymmetric models 



In the previous paragraph we derived a complicated network of diffusion equations (|6.57|) 
that couples all particle species in the hot plasma. In principle, after specification of decay 
rates, diffusion constants and CP-violating sources, it is possible to solve the transport 
equations numerically. However, analytic progress can be made by using conservation 
laws and neglecting interactions that are slow compared to the relevant time scale. An 
interaction with rate T can be neglected if the typical interaction time is large compared 
to the average time a particle spends diffusing in front of the wall before being caught by 
the wall, which is equivalent to |k| 

D 



2 <r~ 1 . (6.58) 



The electroweak ( "weak" ) sphaleron interaction with rate T ws is slow in precisely the 
sense of ( p.58| ) (unless the wall velocity is particularly small, i.e. v w <^0.01). In the 
following we will therefore assume baryon and lepton number conservation and include 
the weak sphalerons only at the end of the calculation. The neglect of the weak sphalerons 
allows us to completely forget about leptons in our transport equations and compute only 
the quark and Higgs densities. 

The processes we do take into account are the supergauge interactions, the strong 
sphaleron interactions, and those described by the Lagrangian 

Ant = ytt c q 3 H 2 + ytt c q 3 h 2 + y t t%h 2 - y t ^i c *q^H x + y t A t t c q 3 H 2 + h.c. (6.59) 

Via terms of type \sh\H 2) the singlino, s, is coupled to the quark-Higgs system. In the 
case of transitional CP- violation triggered in the singlet sector (see fig. [?p, the singlino 
receives the largest CP-violating source term of all particles, even though its impact on 
baryogenesis is suppressed by the small coupling A. For the sake of simplicity we ignore 
this interesting contribution in the following. We assume the supergauge interactions to 
be in equilibrium. The chemical potential of any particle is then equal to that of its 
superpartner, with exception of the singlet field. It is convenient to define the chemical 
potentials \i v = (// u c+// fic )/2, /i Ql = (/4t+//d+//fi+/^)/4, /j, Hl = (fi H o+[i H -+[i~ h o+fi~ h -)/4:, 
etc. In this notation the interaction terms take the formQ 

(T y + Y yA )(^ H2 + /iQ 3 + /x T ), T yil (n Hl - fi Qs - /i T ), 

T ss (2jj,Q 3 + 2/j,q 2 + 2/i Ql + fi T + /i B + fi c + /is + + I^d), 

IV (jJ^Hi + Hh 2 ) , r m (jmq 3 +fj, T ), r Hl fi Hl , T H2 fi H2 . (6.60) 

13 In contrast to refs. |38],|64|] we count all left-handed particles and the corresponding superpartners 
with positive signs. 
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The rates in the first line are related to the interactions ( |6.59| ). T ss denotes the strong 
sphaleron rate. T h f is due to Higgsino helicity flips induced by the fihih 2 term. T Hl2 
and T m correspond to Higgs and axial top number violating processes, present only in the 
phase boundary and the broken phase. 

Because of the small Yukawa couplings of the first and second family quarks, these 
particles are in very good approximation only produced by strong sphalerons. Hence their 
number densities are algebraically constrained. If the system is near thermal equilibrium, 
number densities and chemical potentials are related by 

rii = \kiiiiT 2 (6.61) 
6 

where ki is the appropriate sum over statistical factors n introduced in ( |6.54| ), e.g. k,Q l = 
N c (k u + K d + Ku + kj), ku = N c (k u c + Ky,c), k Hl = (k h q + k h - + K~ h o + n~ h -) , etc. N c = 3 



denotes the number of colors. In the massless limit used in ref. |)4j one obtains kq ± 2 3 = 1<: 
ku = ko = ■■■ = kx = 9, ks 12 = 6. Using baryon number conservation the strong 
sphaleron rate reads 



r s S (2/iQ 3 H h Hd) = r ss 



2 + 9^)^ + (l- 9 |)^ 



(6.62) 



To arrive at this expression we assumed that all the squark partners of the light quarks 
are degenerate in mass. Assuming equilibrium for the strong sphalerons we obtain 

2k B + 9k Q3 

» T = 9k T -k B ^" (6 - 63) 

The validity of this assumption will be discussed below. 

We are now able to write down the reduced set of diffusion equations for the relevant 
particle species Q 3 , Hi and H 2 

- AV q n Q . A + (T y + T yA )[fi H2 + BfiQ.j] - T yfl [fi Hl - BfiQ 3 ] + BT m fi Q:j = S Qa 
—kn^h^H! + Fy^lfiH! — Bfig 3 } + T h f[fi Hl + fi H2 ] + T Hl fi Hl = Sjji 
-k H2 V h ^H 2 + (T y + Y yA )[^H 2 + B/j,q 3 ] + T hf [j2 Hl + fiH 2 ] + T H2 j2 H2 = Sh 2 (6.64) 



where 



.4 



9k T k Q3 + 9k B k Qs + 4k B k T 
9kx — k B 



B = kB + 9kT + 9kQ > (6.65) 
9k T — k B K J 

and T>i = Di-^ +^^- These equations result from summing the diffusion equations of 
particles which belong to the same color and SU(2) multiplets. We have taken a common 
diffusion constant, D q , for the (s)quarks, as well as one for the two Higgs doublets, Dh- 
Notice that the effects of hypercharge screening have been neglected, which can be shown 



to affect the created baryon asymmetry at most by a factor of order one . 

We keep the rates related to the top Yukawa interactions finite. If these interaction 
are in equilibrium, the resulting diffusion equations are sourced only by the combination 
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Sh ± — Sh 2 , because of the constraint j2 Hl +^h 2 = 0- As a result the dominant contribution 
to the chargino source terms cancels, because the corresponding terms for K[ and h\ are 
exactly of the same size. This would not be true for the 7', 5' contribution if we used the 
dispersion relations for canonical momenta. 

In the MSSM the full diffusion equations (|6.64|) have already been studied in ref. [[FT 



and later in ref. [73|. Applying the closed time path formalism to calculate the source 



terms, the diffusion equations (|6.64|) were analyzed recently in ref. 174 



6.4 The baryon asymmetry 

6.4.1 Solution of the diffusion equations 

In the previous section we derived differential equations for the chemical potentials of 
the various particle species contained in the hot electroweak plasma. Baryon number 
violation has been neglected throughout this calculation. However, what we set out 
to compute was the total baryon asymmetry created during the phase transition. So 
before solving the network of diffusion equations, we turn to baryon number generation 
by weak sphaleron processes, which is fueled by the chemical potential of left-handed 
quarks, /ig L = jiQ 1 + /ig 2 + hq z . Using baryon number conservation and eq. ( |6.63|) we find 



fc Q 3 + 2k T 
9k T - k 



B 



2k B 2k B 



k, 



Ay 



fiQ 3 = CfiQ 3 . 



The evolution of the baryon number density n B is governed by 

-V q n B + 3Q(z)T ws (T 2 u jBl - an B ) = 0, 



(6.66) 



(6.67) 



where we have assumed identical diffusion constants for all quarks and squarks, and 
neglected contributions of leptons. The position dependence of the weak sphaleron rate is 
modeled by a step function Q(z): anomalous baryon number violation is unsuppressed in 
the symmetric phase (z > 0) and suddenly switched off in the broken phase {z < 0). The 
second term in Eq. ( |6.67| ) describes damping of the baryon asymmetry by weak sphalerons 
in the symmetric phase. The parameter a depends on the degrees of freedom present in 
the hot plasma. Taking only the right-handed stop to be light gives a = 48/7 
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From Eq. |6.67| one can easily obtain the baryon to entropy ratio in the broken phase 

135r 



Vb = — 



s 



dz/j 



B T \z e 



u 



(6.68) 

) E3 



27r 2 g*v w T 

where we have taken the entropy density s = (2ir 2 g 3f /45)T 3 and v 

~ 126 is the effective number of degrees of freedom at the phase transition temperature. 
Eq. (|6.68| ) shows that the integral over the left-handed quark number, ul oc \i Bh in the 
symmetric phase determines the final baryon asymmetry. 



We now return to eqs. (|6.64|) in order to compute [i Bh . These linear second order 
differential equations can be solved by finding the appropriate Green's function. We keep 
the discussion general and consider the following set of iV coupled diffusion equations 



/ -k u V n + T U ■■■ -k 1N T> 1N + T 1N \ 



\ -k N1 v N1 + r 



iVl 



—k NN V NN + T NN J 



I Pi \ 



\ J 



( Si \ 



\S N J 



(6.69) 
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where V ab = D ab -^ + v wj=- The corresponding boundary conditions read /i a (M — > oo) = 
0. The matrix valued Green's function G ab is defined by ^2^ =1 (—k ac V ac + T ac )G cb (z) = 
8ab$(z)- in the transport equations (|6.64| ) also position dependent rates are present, 
e.g. T m . They typically vanish in the symmetric phase and become maximal in the 
broken phase. In order to keep the problem analytically tractable we simply model the 
position dependence of these rates by step functions, i.e. T ab {z) = T +a bQ(z) +r„ a ;,6(— z). 
The general structure of the Green's function then reads 



Gah(z) 



©(*)E 



TV JV 

c i+ab e- x * +z + G(-z) (k-a b e~ Xi - M . 

i=l »=1 



(6.70) 



The constants \i± can be computed from det[k ab (— D ab \\ +v w X±) +r ±a6 ] = 0. Of course, 
this procedure supplies us with 4N solutions. However, 2N of them (A$+ < for z > and 
Aj_ > for z > 0) correspond to exponentially growing solutions of eq. ( |6.69 ) and have 
to be discarded. The coefficients Ci± ab can then be computed from the 2N 3 dimensional 
set of linear equations 



N 



^ y [ - k ab D ab \ 2 i+ + k Aj+ + T +ab I Ci +b . 
6=1 

N 

N 

^ \pi+ab Q— ab\ 



6=1 



i—bc 



i=l 



N 



^ ] k a bD ab [Aj_|_Cj + b c — Aj„Cj_;, c ] 
i,6=l 



(6.71) 
(6.72) 
(6.73) 
(6.74) 



where T± ab denote the rates for positive and negative z. Eqs. ( |6.71| ) and ( |6.72| ) result 
from solving the homogeneous version of eq. (|6.69|) in the range of positive and negative 
z, respectively. Only 2N 2 (N — 1) of them are independent. The continuity of the Green's 
function at z = is guaranteed by eq. ( |6.73|) . Finally, eq. ( 6.74Q is obtained by integrating 
the definition equation for the Green's function on an infinitesimally small interval around 
z = 0. Diffusion constants D ab and statistical factors k ab differing in broken and symmetric 
phase can be treated along the same lines. 

Once the Green's function is known one can easily compute the chemical potentials. 
Applying the general formulas to the diffusion equations ( |6.64| ) we identify H\ = fj,Q 3 , 

= Mifn A*3 = A*ir 2 > Si = S Qa , S 2 = S Hl and S 3 = Sh 2 - Using eqs. ( |6.68| ) and ( |6.66| ) we 
obtain for the baryon to entropy ratio 



i35r, 



2n 2 g*v w T 



dze~ vl 



dz'Gi a (z-z')S a (z'). 



(6.75) 



Since the Green's function consists only of exponentials, the z-integration can be per- 
formed analytically. The evaluation of this expressions is performed in the next paragraph. 
Also the various approximations we used in its derivation will be discussed. 



34 



6.4.2 Numerical evaluation and discussion 



Before starting to calculate a numerical value for the emerging baryon asymmetry we 
discuss the validity of the assumptions and approximations made in the derivation of 
eq. ( |6.75D . This requires the specification of the various parameters that enter the diffusion 



equations. For the diffusion constants of quark and Higgs fields we take ]73] 



D 9 =~, D h =™. (6.76) 



We use the rates [[T4|] 



T y + Y yA = 0.015T, r yM = 0, T m = 0.05T6(-^) 

Y hf = 0.016T, T Hl = T H2 = 0.0036T6(-z). (6.77) 



T yfi is strongly suppressed because it involves heavy left-handed stop states ( |6.59 ). T 



Th 1 and Tjj 2 are present only broken phase. We model their position dependence by a 



step function. The weak and strong sphaleron rates are given by J76| 



T ws « 6a*T « 2.2 x 1(T J T 

T ss « 1500r ws « 0.033T. (6.78) 



It has been shown recently that parametrically T^s = C ln(l/g 2 )a4 |ff7| , but lattice 
measurements of the rate are consistent with C ~ l/oc w . The thickness of the bubble 
wall varies considerably in the NMSSM, 1/T < L w < 20/T, i.e. the wall may become 



much thinner than in case of the MSSM, where 20/T < L w < 30/T has been found |35]. 
Calculations of the wall velocity in the SM lead to 0.36 <v w < 0.44 [|67|]. Gauge fields in 
the hot plasma diminish this result . In supersymmetric models there arise additional 
friction terms from the SUSY particles, in the first place from a light top squark. In the 
MSSM, this can bring down the wall velocity to v w ~ 10~ 2 f68f . In the following we treat 
v w as a free parameter and examine its impact on the emerging baryon asymmetry. 

Let us now summarize the approximations leading to eq. ( |6.75| ) for the baryon to 
entropy ratio: 

• Assumption 1: L w > ^ so that most particles in the plasma are indeed accurately 
described by the WKB approximation in their interaction with the bubble wall. For 
typical wall thicknesses in the NMSSM, 5/T < L w < 15/T, this assumption is very well 
justified, although for the very thinnest walls L w ~ 1/T the WKB approximation becomes 
questionable. 

• Assumption 2: v w < j^, the thermalization condition that guarantees the appli- 
cability of the fluid ansatz. In case of (s) quarks this condition is satisfied for v < 0.4 
even for rather thin walls with L w ~ 5/T. Higgs particles thermalize much slower, as 
can be deduced from their large diffusion constant. Even for the largest wall thicknesses, 
L w ~ 20/T, rather small velocities, v w < 0.1, are required. 

• Assumption 3: v w < We work to linear order in the wall velocity, which is 
only justified if the wall moves slower than the speed of sound in the plasma. Otherwise, 
diffusion in the region in front of the wall, giving rise to "non-local" baryogenesis, is no 
longer possible. However, in case of v w > -j= the fluid approximation would break down 
anyway, i.e. from assumption 3 no new constraints result. 
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• Assumption 4: r„, s < v^/D, so that the back- reaction of the baryon number violating 
processes can be neglected in the diffusion equations Q6.64] ). According to eqs. ( |6.76D and 
( |6.78| ) this requires v w > 0.01 (quarks) and v w > 0.04 (Higgs particles). 

• Assumption 5: r ss > v 2 w jD in order to put the strong sphaleron interaction to 
equilibrium. In this case the strongest constraint, v w < 0.45, result from the quarks, 
which is easily satisfied. 

The calculation of the final baryon asymmetry still requires the specification of the 
source terms (|6.57| ) which enter eq. ( |6.75|) . We concentrate on the source terms induced 



by the semi-classical force and neglect those related to spontaneous baryogenesis. We 
include the source terms of the top quark ( |6.32j ) and the charged Higgsinos ( |6.40|) . For 



sake of simplicity we neglected the neutral Higgsinos, which should affect the final baryon 
asymmetry at most by a factor of order one. Also the singlino is disregarded. 

In order to evaluate the source terms ( |6.57 ) the bubble wall profile is required. We 



approximate our numerical solutions by a kink-ansatz with a common wall thickness for 
all fields present in the bubble. However, we allow for a variation in tan {3 in the bubble 
wall by taking \H%(z)\ = \H%(z)\ (tan fa + const • \H^(z)\ 2 \), in the spirit of fig. |a. The 
relevant source terms are obtained from combining eqs. ( |6.32|) , (|6.40 ) and fl6.57f ) 



St = r^{(A)9>;+(B){rnl)>0>)i (6.79) 

s hl = 7 |^(^)(^-7b in2 K) + ^^ 2 (M) , + 

+ (B){m\y(e' }i - 7 ^ in2 K) + ^ 2 (h))) ■ (6.so) 

S~ h2 = S~ hi (6.81) 



with 



Here iV c = 3 is the number of colors and E 2 = p 2 + m?, E 2 = p 2 + m 2 and p\ = p 2 — p 2 . 
With exception of {p 2 z /E) all thermal averages in eqs. ( |6.79| ) - ( |6.81| ) are performed with 



massive distribution functions for fermions ( |6.54|) , which ensures the decoupling of heavy 



particles. Taking into account only the source terms of eqs. ( |6.79| ) - ( |6.80| ), we have 
Sq 3 = S t , S Hl = S^ and S H2 = S^. 

In the approach with canonical momenta (which we had in the first version of this 
paper), there was a source Sj for the stop and S^ ^ S^ 2 . The latter leads to an important 
contribution to the Hi — if 2 -combination, the well-known term proportional (tan 0)' and 
to an additional term for the singlet field. These are absent now. The brackets (A) and 
(B) in the approach with canonical momenta are 

(An = (^f^) + , (B-> = {^^) + (6-83) 

For p^> m they differ by a factor from the ones for kinetic momenta ( |6.82|) . 
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The generated baryon asymmetry depends on the squark spectrum because of the 
potential suppression due to strong sphalerons. In the following we consider four different 
squark spectra which are listed in table |[ Case A corresponds to the massless limit, 
where strong sphaleron suppression is most efficient. In case B all squarks are assumed to 
be heavy, with exception of the right-handed stop, which is taken massless. The spectra 
C and D are obtained from the universal pattern of SUSY breaking discussed in section 
2. The first and second generation squarks are almost degenerate in mass with b 2 , thus 
2k B = kq 1 = kq 2 = 2ku = . . . Case C corresponds to the example of explicit CP-violation 
discussed in context of fig. [| while case D is the squark spectrum of the transitionally 
CP-violating example of fig. [7]. 

According to eq. (|6.75 ) the produced baryon asymmetry is proportional to the param- 



eter C, which is related to the chemical potential of left-handed quarks ( |6.66|) and also 



given in table ^. It turns out that the results for different spectra can indeed be obtained 
by rescaling with the relevant C parameters, i.e. the indirect impact of the statistical 
factors in the diffusion equations ( |6.64| ) is small. For the massless case A, and more 
general for the case of degenerate squarks, C vanishes. Baryon production is completely 
suppressed by rapid strong sphalerons transitions, which is a well known result |75| . The 



cancellation disappears for non-degenerate squarks. Even for the realistic spectra C and 
D, resulting from universal SUSY breaking, there is only a mild suppression by a factor 
of 3-4 relative to the idealized spectrum B (usually assumed in context of the MSSM). 

To begin with let us consider the source term of the top quarks (|6.79| ), which in the first 
place depends on the variation of the phase of the top quark mass along the wall, S t oc 59 1 , 
and the wall thickness L w . From eqs. ( |5.14 ) and ( |5.22| ) we obtain 58 t = 502 = cos 2 (j3)58, 



where 62 denotes the phase of and 9 is the (gauge invariant) sum of the phases of 
the two Higgs fields, that entered the bubble equations in section 5. As a consequence, 
St is rather sensitive to the Higgs vev ratio. Baryon production from the top quark 
source is particularly efficient for small values of tan(/?T). In fig. |9| we present the baryon 
asymmetry induced by S t as a function of the wall thickness and velocity, where we used 
v c /T c = 1.60, T c = 101 GeV and y t = 1.015 (corresponding to tan/3 = 5). 

Our results^ would be slightly enhanced by larger values of v c /T c . We measure the 



14t 



In numerical evaluations we still used the brackets (A can ), (B can ) of eq. faM) instead of (A), (B) of 
eqs. ( 6.82| ) in order not to be forced to completely repeat the whole analysis. This introduces a deviation 



of order (1), in the case \p\ ^> m a factor ^ instead of ^4 in the brackets (A) and (B). This is well 
within the accuracy of the present investigation. 
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0.385 
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125 


-100 


418 


238 


334 


375 


7.26 


4.56 


3.60 


7.20 


0.125 


D 


100 


-100 


371 


201 


217 


325 


8.43 


4.98 


3.87 


7.74 


0.102 



Table 2: Squark spectra and corresponding statistical factors used in the discussion of the 
baryon asymmetry. For the sets C and D we used A = 0.05, k = 0.4, m = 200 GeV and 
tan/? = —5. Furthermore, we took x = —100 (—150) GeV in case C (D). The sets A and 
B do not follow the universal pattern of SUSY breaking considered in section 2. (Units 
in GeV.) 
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Figure 9: Top quark contribution to the baryon asymmetry for the spectrum D in units 
of 2 x 10 -11 normalized by 59t as a function of the wall velocity for different values of the 
wall thickness L w = 20/T, 10/T, 5/T, 3/T (from below). 



baryon asymmetry in units of 2 x 10 , which is the lower observational bound |69| 



"* = 2x~l(F"' (6 ' 84) 



In fig. ^ we display the dependence of the baryon asymmetry on the wall velocity for 
different values of the wall thickness L w = 20/T, 10/T, 5/T, 3/T. We used the squark 
spectrum D in the evaluation. Results for other spectra can be easily obtained by rescaling 
with the relevant C parameters of table |^. For L w <r ^10/T baryon production is most 
efficient for v w ~ 0.02 — 0.03. Interesting enough, this is precisely the range of wall 
velocities found in recent calculations in the MSSM [^,|6^]. rj becomes small for large 
values of v w because transport in front of the wall gets more and more inefficient. For 
small wall velocities deviation from equilibrium becomes too small to generate a relevant 
baryon asymmetry. We observe an approximate dependence in the generated baryon 
number. This behavior is expected from the explicit expression for the source term ( |6.79|) , 
which contains three derivatives with respect to z. Since the calculation of rj ( |6.75| ) requires 
one (numerical) integration of the source over z', approximating <9 g ~ 1/L W we obtain 
7] oc 1/T^. As a result baryon production from the top quark source is most efficient 
in the case of slow (v w ~ 10~ 2 ) and thin walls, and for non-degenerate squark spectra. 
It turns out that the Higgs source term behaves in a similar way. In fig. |9| we assumed 
5/3 = 10~ 3 in the bubble wall which is a typical value according to section 4. However, 
even large values, 5(3 = 10~ 2 , change the result only by a few percent. From fig. ^ we 
approximately obtain for the maximal value of the baryon asymmetry 

v ^^ A7 {i~t) 6dt - (6 - 85) 

However, for medium and large Higgs vev ratios even small values of 56 t ~ 0(1O~ 2 ) are 
difficult to achieve, since 58 t = cos 2 (/3)56. In the case of tan/5 = 5, for instance, we obtain 
56 t = 56/26, which leads to a significant suppression. In the example of transitional 
CP- violation (fig. 0) we have tan/3 = 5, L w = 3 and 59 ~ 1/20 leading to ^ ax ~ 1/4. 



38 




t 

V* 



- 




A (b) - 

/ \ 
/ \ 






/ \ : 
/ ' N \ 






/ / \ \ 






/ / \ \ 






// \ \ 





















log 



10 v 



Figure 10: The chargino contribution to the baryon asymmetry in units of 2 x 10 -11 
as a function of the wall velocity for different values of the wall thickness L w = 
20/T, 10/T, 5/T, 3/T (from below). We use the squark spectrum C and the example 
of explicit CP-violation considered in the context of fig. |||. (b): The same quantity for 
the transitionally CP- violating bubble wall of fig. [7] and the squark spectrum D. 



In case of explicit CP-violation, the top quark contribution to the baryon asymmetry is 
much smaller even, and can be safely neglected. 

We now come to the charged Higgsino contribution to the baryon asymmetry fu- 
eled by the source terms ( |6.80| ) and (|6.81| ) for the examples of explicit and transitional 



CP-violation already discussed in the context of figs. |6] and |7|. Again, we present our 
results 14 as a function of the wall velocity for different values of the wall thickness 
L w = 20/T, 10/T, 5/T, 3/T. Concerning the dependence of the generated baryon asym- 
metry on the wall velocity and the wall thickness we find a similar behavior as in the case 
of the top quark source: rj becomes large for thin walls and v w ~ 10~ 2 . 

The baryon asymmetry generated from the chargino dispersion relation is shown in 
fig. [l(| If the top Yukawa interactions were in equilibrium, as was assumed in elder work 
on the subject, this contribution is completely erased, because of the equal source terms 
for both Higgsinos ( |6.80|J6.8T| ). Even with finite top Yukawa rates no baryon asymmetry is 



generated from this source, if both Higgs fields have equal rates in the diffusion equations 
( |6.64| ), i.e. if T y + T v a = Y yfl and T Hl = T H2 . In our scenario an asymmetry in the top 
Yukawa rates is inevitably induced by the heavy left-handed stop, leading to a strong 
suppression of T yfi ( |6.77| ). In the evaluations of fig. [L(] we used 8f3 = 10~ 3 . If we increase 
the change in the Higgs vev ratio by a factor of 10, the results only change by about 15 
percent. 



From fig. 10 a for the case of explicit CP- violation we read of 



^-~°- 8 (t^) sin ^)- (6-86) 

Thus CP- violating phases of order 10 _1 are required to account for the observed baryon 
asymmetry in this scenario. Phases of this size are compatible with the EDM experiments 
only if the first and second generation squarks are heavy or in the case of accidental 
cancellations. Our findings for baryogenesis from explicit CP-violation in the NMSSM 



39 




0.75 



0.25 



75 125 175 225 275 325 375 




150 200 250 300 350 400 450 



M 2 [GeV] -> m Ua [GeV] -> 

Figure 11: Chargino contribution to the baryon asymmetry in units of 2 x 1CT 11 as a 
function of the gaugino mass parameter M 2 (a) and as a function of the soft scalar mass 
parameter mjj 3 . In both cases we otherwise used the parameter set of explicit CP- violation 
(fig. |), </v = 0.1, L w = 5/T, v w = 0.03 and 5(3 = 10~ 3 . 



are very similar to the results in the MSSM 73 . Still, baryogenesis in the NMSSM is 



slightly more effective because of the thinner bubble walls. 

The scenario of transitional CP-violation is a very attractive alternative, being not 
in conflict with the EDM experiments. As shown in fig. |10|a, it can generate the baryon 
asymmetry for a considerable range of wall velocities, especially for L w = 3/T. It even 
allows for a certain baryon "over-production" to compensate for a partial washout of the 
generated baryon asymmetry due to the existence of bubbles with opposite CP-properties. 
The introduction of small explicit phases lifts the energetic degeneracy of the two sort 
of bubbles, leading to a more rapid nucleation of bubbles of the true vacuum and a net 
baryon production. 

Up to now we centered the discussion of baryon production around two particular 
sets of NMSSM parameters. However, the results may considerably change if different 
regions in the parameter space are considered. Of course, baryon production becomes 
inefficient if the particles which supply the CP-violating source terms, in the first place 
the Higgsinos, are heavy and decouple from the plasma. The thermal averages, which 
enter the sources ( |6.79|) - Q6.81]) , already decrease by a factor of 10 if m/T is raised from 
one to five. In case of m/T = 10 the suppression is of order O(10~ 3 ) and successful 
baryogenesis becomes extremely difficult. If universal SUSY breaking is assumed this 
occurs in the regime M > 700 GeV. 

The chargino source term is rather sensitive to the relative size of the Higgsino mass 
parameter, fi, and the mass of the SU(2)-gauginos, M 2 . In fig. |TT| a we present the depen- 
dence of the generated baryon asymmetry on M 2 , while the remaining parameters as well 
as the bubble wall profile are taken from the example of explicit CP-violation considered 
in fig. [|, i.e. the effect of a change in M 2 on the bubble wall is ignored. Moreover, we use 
the squark spectrum C, 4>n = 0.1 L w = 5/T and v w = 0.03. We find a resonance structure 
for r], with the peak located at M 2 ~ = 226 GeV. Fig. [10| a corresponds to M 2 = 103 
GeV. In fig. [Tl| b we present the dependence of the baryon asymmetry on the soft breaking 
scalar mass parameter mu 3 - Even though this parameter has no direct impact on the 
chargino source terms (|6.80|) , ( |6.81|) it alters the statistical factors kQ 3 and kx due to 
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varying stop masses. Baryogenesis becomes less efficient if the two stop stated are more 
or less degenerate. In the case under consideration this happens for mu 3 ~ rriQ a = 329 
GeV. In fig. [T0| a we used my. A = 256 GeV. 



7 Conclusions 

In our computation we used the semi-classical approximation to describe the interaction of 
the particles in the plasma and the propagating bubble wall which in context of the MSSM 
has been introduced in ref. fl64| . Our work differs in various aspects from that analysis. 
We included contributions from top quarks which have previously been neglected and 
took into account also position-dependent mixing in the chargino mass matrix ("flavor 
contribution"). In ref. different results were obtained due to an error in the transport 



equations, which prevented the cancellation between the dominating helicity parts in the 
and h\ source terms. With dispersion relations for canonical momenta, the charged 
Higgsinos provide in the MSSM CP-violating source terms which are proportional to the 
variation of tan/? in the bubble wall. This result would agree well with refs. p9"l|5^, [58)1 
where different methods where used to determine the CP- violating source terms. Recently 
5 f3 independent contributions to the source terms in the MSSM have been calculated ||60|| . 
However, these arise from higher orders in a Higgs insertion expansion and should not 
be included in our semi-classical approach. In the NMSSM non-vanishing source terms 
for stops and charginos (and for the top quark) are generated even for constant tan/3 
because of the changing singlet field, and due to CP-violating bubble walls. These arise 
either from explicit CP-violation, which in case of the NMSSM is possible already in the 
tree-level Higgs potential, or from transitional CP-violation. 

Using kinetic momenta in the dispersion relations as suggested convincingly in ref. |j4 



hi and h 2 sources are exactly equal and tan /3' effects cancel. Thus one has to consider a 
source symmetric in the Higgses. This contribution can be also sizeable, if the Yukawa 
interactions are not in equilibrium |6^] and if the left-handed stop is heavy [[^,[74]]. This 



statement ist unchanged in our singlet model. The WKB results and their difference to 
other work should be further discussed using quantum transport equations. 

In our work we investigated the impact of different squark spectra on the emerging 
baryon asymmetry, while in ref. pi only the massless approximation was considered. As 



a consequence, we obtained a non-vanishing baryon asymmetry, while assuming the strong 
sphalerons to be in equilibrium. Baryogenesis becomes most efficient for non-degenerate 
squarks. The generated baryon asymmetry is roughly proportional to Since in the 

NMSSM the bubble walls can be considerably thinner than in the MSSM, more baryons 
can be produced. The baryon asymmetry strongly depends on the velocity of the bubble 
wall, which in the NMSSM has not been computed so far. Interesting enough, we find 
that the maximal baryon number is generated for v w ~ 10~ 2 . This is a typical value for 
the wall velocity in the MSSM. 

Finally, we studied the dependence of the chargino source terms on /x, M 2 ,mg 3 and 
mfj . We found an enhancement if \x and M 2 become degenerate. This resonance be- 



havior has also been reported in refs. |4^,[52|,|6^,|64[]. If uiq ~ m j/,, baryon production 



is partially suppressed by strong sphalerons. Concluding, the NMSSM can account for 
the observed baryon asymmetry, especially if \x ~ M 2 and if the squark spectrum is 
rather non-degenerate, or if transitional CP-violation occurs. Both scenarios point to a 
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non-universal pattern of SUSY breaking. 
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